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1  Fully  Developed  Laminar  Flow  in  Trapezoidal  Grooves 
1.1  Introduction 

Internally  grooved  ducts  are  used  in  process  equipment  to  improve  heat  transfer 
during  condensation  and  evaporation.  In  some  cases,  such  as  a  refrigeration  cycle 
evaporator  or  condenser  with  internal  grooves,  the  vapor  flow  is  cocurrent  with  respect 
to  the  liquid  flow.  In  axially  grooved  heat  pipes,  the  vapor  flow  is  countercurrent  to  the 
liquid  flow.  The  interfacial  shear  stress  due  to  the  cocurrent  or  countercurrent  vapor 
flow  contributes  to  the  liquid  pressure  drop,  which  can  significantly  affect  the  heat 
transfer  capacity  of  the  grooved  surface.  The  objective  of  the  present  research  was  to 
numerically  model  the  flow  of  liquid  in  trapezoidal  grooves  using  a  finite  difference 
approach  in  order  to  provide  accurate  information  on  the  effects  of  cocurrent  and 
countercurrent  vapor  flow  on  the  pressure  drop  in  the  liquid.  This  geometry  was 
chosen  due  to  the  fact  that  rectangular  and  triangular  grooves  are  special  cases  of  the 
trapezoidal  groove,  thus  making  the  analysis  as  general  as  possible. 

DiCola  (1968)  solved  the  conservation  of  momentum  equation  for  the  laminar  flow 
of  liquid  in  rectangular  grooves  with  a  uniform  shear  stress  imposed  at  the  liquid- 
vapor  interface  using  separation  of  variables.  Unfortunately,  this  manuscript  is  no 
longer  available  in  the  open  literature.  Schneider  and  DeVos  (1980)  provided  the 
exact  solution  determined  by  DiCola  (1968),  along  with  an  expression  for  the  friction 
factor  which  approximates  the  exact  solution  to  within  1%  by  using  the  first  term 
of  the  infinite  series  solution.  This  expression  was  used  by  Schneider  and  DeVos  to 
determine  the  nondimensional  heat  transport  capacity  of  axially-grooved  heat  pipes. 
Upon  examination  of  the  DiCola  equation,  it  is  obvious  that  the  rectangular  groove 
is  completely  full;  i.e.,  the  meniscus  contact  angle  is  (j>  —  90°. 

Ayyaswamy  et  al.  (1974)  used  the  Galerkin  boundary  method  to  solve  the  Poisson 
equation  to  determine  the  fluid  velocity  in  triangular  grooves.  In  this  study,  interfacial 
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shear  stress  was  zero,  and  the  groove  half-angle  and  contact  angle  were  varied  from 
9  =  5  to  60°  and  <fi  =  0.1°  to  the  full  groove  condition  (6  +  0  =  90°).  The  results 
were  presented  graphically  and  in  tabular  form,  which  included  the  cross-sectional 
area,  mean  velocity,  hydraulic  diameter,  and  Poiseuille  number.  It  was  found  that 
the  Poiseuille  number  increased  monotonically  with  meniscus  contact  angle. 

Ma  et  al.  (1994)  determined  the  Poiseuille  number  for  the  flow  of  liquid  in  trian¬ 
gular  grooves  with  liquid-vapor  frictional  interaction.  The  groove  half-angle  ranged 
from  20  <  9  <  60°  and  the  meniscus  contact  angle  was  varied  from  (f>  =  0  to  60°. 
The  conservation  of  momentum  equation  was  transformed  into  the  Laplace  equation 
in  terms  of  liquid  velocity,  which  was  solved  using  separation  of  variables  and  lin¬ 
ear  superposition.  Difficulties  were  encountered  with  respect  to  application  of  the 
liquid-vapor  interface  boundary  condition,  since  the  liquid-vapor  interface  velocity 
is  an  unknown  function  of  the  vapor  velocity.  The  methodology  used  to  overcome 
this  difficulty  required  an  experimentally  determined  coefficient.  A  dimensionless 
liquid- vapor  interface  flow  number  was  introduced  to  account  for  the  interfacial  shear 
stress.  This  value  determined  the  relative  velocities  of  the  liquid  and  vapor  at  the 
liquid-vapor  interface,  which  impacted  the  magnitude  of  the  Poiseuille  number.  It 
was  found  that  the  friction  factor  increased  with  the  interface  flow  number  and  con¬ 
tact  angle.  Results  from  the  experiment  executed  by  Ayyaswamy  et  al.  (1974)  for  no 
liquid- vapor  shear  stress  showed  an  excellent  comparison  with  the  analytical  solution 
over  the  ranges  of  groove  half-angle  and  contact  angle  mentioned  above. 

Romero  and  Yost  (1996)  analyzed  the  flow  of  liquid  in  a  triangular  groove  with  no 
shear  stress  at  the  liquid-vapor  interface.  In  particular,  the  flow  from  a  sessile  drop 
into  the  groove  was  of  interest.  A  nonlinear  partial  differential  diffusion  equation  was 
presented  which  described  the  time-dependent  height  of  liquid  in  the  groove  in  terms 
of  the  groove  geometry,  meniscus  contact  angle,  and  fluid  properties.  A  simplified 
similarity  solution  was  presented  for  the  region  which  was  far  from  the  sessile  drop.  A 
full  similarity  solution  was  also  shown  which  accounted  for  the  region  near  the  fluid 
droplet.  It  was  found  that  the  wetting  front  position  was  proportional  to  (Dt)1/2, 
where  the  diffusion  coefficient  D  was  related  to  the  groove  geometry,  liquid  viscosity, 
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and  liquid- vapor  surface  tension. 

Lin  and  Faghri  (1997)  modeled  the  flow  of  liquid  in  the  triangular  grooves  of 
a  rotating  miniature  heat  pipe.  A  correlation  for  the  friction  factor  was  provided 
in  terms  of  the  shear  stress  at  the  liquid-vapor  interface.  The  laminar  flow  in  the 
triangular  groove  was  solved  using  a  finite  element  technique  for  side  lengths  ranging 
from  hV  1  +  tan2  9  =  0.2  to  0.65  mm  and  liquid-vapor  shear  stress  riv  =  7.7  x  10-5  to 
0.055  N/m2  for  a  groove  half-angle  of  9  -  20°  and  meniscus  contact  angle  <j>  =  30°. 
A  regression  analysis  was  used  to  represent  the  data  to  within  ±2.8%. 

Khrustalev  and  Faghri  (1999)  analyzed  the  fully  developed  laminar  flow  of  liquid 
and  vapor  in  miniature  heat  pipes  using  a  finite  element  solution.  In  particular,  the 
case  in  which  the  vapor  velocity  was  high  and  the  cross-sectional  areas  of  the  vapor 
and  liquid  were  comparable  was  of  interest.  It  was  assumed  that,  with  respect  to 
the  vapor  flow,  the  liquid  velocity  at  the  liquid-vapor  interface  was  zero.  For  the 
liquid  flow,  the  shear  stress  at  the  interface  was  equal  and  opposite  to  that  of  the 
vapor.  This  meant  that  the  velocity  gradient  in  the  liquid  was  related  to  that  of 
the  vapor  via  a  ratio  of  absolute  viscosities.  The  momentum  equation  for  the  vapor 
domain  was  first  solved  to  determine  the  shear  stress  distribution  at  the  liquid- vapor 
interface.  Then  the  momentum  equation  for  the  liquid  domain  was  solved  using  the 
shear  stress  information  from  the  vapor  solution.  The  results  presented  were  for  a 
specific  heat  pipe  geometry  that  matched  a  previous  experimental  study.  It  was  found 
that  the  shear  stress  at  the  liquid- vapor  interface  was  not  uniform,  being  greater  near 
the  point  of  contact  with  the  solid  groove  wall.  This  effect  was  more  significant  for 
smaller  values  of  meniscus  contact  angle.  In  addition,  the  shear  stress  at  the  interface 
became  more  uniform  as  the  vapor  space  became  more  restricted. 

Kolodziej  et  al.  (1999)  analyzed  the  gravity-driven  flow  of  liquid  in  a  triangular 
groove  with  no  shear  at  the  liquid-vapor  interface.  The  shape  of  the  liquid-vapor 
interface  was  determined  in  terms  of  the  Bond  number  and  meniscus  contact  angle. 
Starting  with  the  Young-Laplace  relation,  a  nonlinear  boundary-value  problem  for 
the  liquid-vapor  interface  shape  was  solved.  The  flow  field  was  then  solved  for  the 
friction  factor  using  the  boundary  collocation  method.  The  range  of  parameters  was 
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as  follows:  groove  half-angle  9  =  5  to  70°,  meniscus  contact  angle  <f>  —  5  to  45  ,  and 
Bond  number  Bo  =  pgb2/c  =  0.001,  0.01,  0.1,  1.0  and  5.0,  where  b  is  the  vertical 
distance  from  the  liquid-vapor  interface  to  the  bottom  of  the  groove.  It  was  found 
that  the  Bond  number  had  a  significant  effect  on  the  friction  factor  of  the  flow. 

The  objective  of  the  present  study  was  to  determine  the  mean  velocity,  volumetric 
flow  rate,  and  Poiseuille  number  for  the  flow  of  liquid  in  trapezoidal  grooves.  The 
effect  of  vapor  flowing  over  the  liquid-vapor  interface  was  accounted  for  by  relating 
the  liquid  velocity  gradient  to  the  friction  factor  of  the  vapor.  This  approach  assumed 
that  the  liquid  velocity  did  not  affect  the  vapor  velocity;  i.e.,  the  vapor  velocity  at  the 
interface  was  zero.  In  addition,  the  variation  of  the  shear  stress  along  the  liquid-vapor 
interface  (Khrustalev  and  Faghri,  1999)  was  neglected.  The  conservation  of  momen¬ 
tum  equation  was  solved  using  Gauss-Seidel  iteration  with  successive  over-relaxation. 
The  analysis  was  validated  using  the  results  of  several  previous  studies  including  the 
flow  of  liquid  in  rectangular  and  triangular  grooves  with  liquid-vapor  interaction.  The 
mean  velocity,  volumetric  flow  rate,  and  Poiseuille  number  are  presented  in  graph¬ 
ical  and  tabular  form  in  terms  of  the  groove  aspect  ratio  (0  <  (3  <  1.5),  groove 
half-angle  (0  <  6  <  60°),  meniscus  contact  angle  (0  <  <j>  <  90°  -  6),  and  dimension¬ 
less  shear  stress  at  the  liquid-vapor  interface  (-0.45  <  r*v  <  5.0).  The  results  were 
used  to  determine  the  effects  of  groove  fill  amount  on  the  capillary  limit  of  revolving 
helically-grooved  heat  pipes.  The  predictions  of  the  improved  capillary  limit  model 
were  compared  to  the  experimental  data  obtained  by  Castle  et  al.  (2000). 

1.2  Mathematical  Model 

A  constant  property  liquid  flows  steadily  in  a  trapezoidal  groove  as  shown  in 
Fig.  1.  A  meniscus,  which  is  assumed  to  be  circular,  comprises  the  liquid- vapor 
interface  (Bo  «  1).  For  fully  developed  laminar  flow,  u  =  w  =  0  and  v  =  v(x,z). 
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where 


The  dimensionless  radius  of  curvature  is  given  by 

((3  +  tan  9) 

/3  sec  <p  for  9  =  0 
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On  the  liquid- vapor  interface,  a  uniform  shear  stress  in  the  y  direction  is  imposed. 

—  =  T,*  :  0  <  x*  <  /3  +  tan  #,  z*  =  (1  +  <f)  -  x/R*2  -  z*2  (8) 

dn* 

The  dimensional  liquid- vapor  interface  shear  stress  can  be  cast  in  terms  of  the  friction 


factor  of  the  vapor. 


Tlv  — 


Pv  {VvY 


fv  for  co  current  flow 


fv  for  countercurrent  flow 


The  Poiseuille  number  of  the  liquid  in  the  groove  is  given  by 


Po  =  /Re 


The  dimensionless  hydraulic  diameter  is 


Dl  = 


2  (3  +  {0  +  tan  #)  (1  +  <f)  -  R*2  cos'1  J  (P  +  sec  9)  for  9  +  $  <  tt/2 
2  (2/3  +  tan#)  (/3  +  sec#)-1  for  #  +  cf>  =  vr/2 


The  mean  velocity  is  defined  as 


=  T*  /  v*dz*dx* 

Al  Jo  Jo 


where  the  dimensionless  cross-sectional  area  is  given  by 

{13  +  (/?  +  tan 6)  (1  +  d')  -  R’2 cos-1  (  4~)  for  6  +  (p  <  tt/2 

\R  )  (13) 

2/3  +  tan  #  for  9  +  <j>  =  7r/2 

See  Appendix  A  for  derivations  of  the  dimensionless  Poisson  equation,  dimensionless 
boundary  conditions,  and  dimensionless  hydraulic  diameter  and  area. 


1.3  Numerical  Model 

The  elliptic  Poisson  equation  given  in  eqn.  (3)  with  mixed  boundary  conditions 
[eqns.  (4),  (5)  and  (8)]  was  solved  using  Gauss-Seidel  iteration  with  successive  over¬ 
relaxation  and  a  second-order  central  differencing  scheme  (Burden  and  Faires,  1985). 
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The  symmetry  boundary  condition  was  modeled  using  a  first-order  forward  difference, 
where  the  shear  stress  boundary  condition  employed  a  first-order  backward  difference. 
The  convergence  criteria  for  the  iterative  solution  was  set  to  e  =  10-8  for  each  case.  A 
grid  independence  check  was  made  in  which  the  number  of  grids  in  each  direction  was 
doubled.  When  the  value  for  the  Poiseuille  number  did  not  change  by  more  than  3%, 
grid  independence  was  considered  to  be  reached.  The  convergence  criteria  was  then 
reduced  by  an  order  of  magnitude  while  maintaining  the  highest  number  of  grids.  If 
the  Poiseuille  number  did  not  change  by  more  than  2%,  the  solution  was  considered 
to  be  independent  of  both  grid  size  and  e.  Otherwise,  a  grid  independence  check 
was  made  at  the  smaller  value  of  e  until  a  converged  solution  was  obtained.  In  fact, 
the  grid  independence  for  423  of  the  446  data  points  reported  was  less  than  1%  as 
documented  in  Appendix  C. 

The  numerical  model  was  tested  against  several  existing  solutions,  such  as  rectan¬ 
gular  ducts  (Fig.  32  and  Table  4) ,  triangular  grooves  without  interfacial  shear  stress 
(Fig.  33  and  Table  5),  and  rectangular  and  triangular  grooves  with  shear  stress  (Figs. 
34  and  35  and  Tables  6  and  7).  See  Appendix  B  for  details. 

Shah  (1975)  determined  the  friction  factors  for  the  laminar  flow  within  ducts  of 
various  cross  sections  using  a  least-squares-matching  technique.  A  comparison  was 
made  of  the  Poiseuille  number  between  the  present  solution  and  those  given  by  Shah 
(1975)  and  Shah  and  London  (1978)  for  laminar  flow  in  a  family  of  rectangular  ducts 
( 6  =  0°,  0.01  <  (3  <  1.0).  The  maximum  difference  was  0.9%. 

Ayyaswamy  et  al.  (1974)  presented  the  friction  factors  obtained  for  laminar  flow 
in  triangular  grooves  using  the  Galerkin  method  of  solution.  Romero  and  Yost  (1996) 
derived  an  equation  for  the  dimensionless  volumetric  flow  rate  of  liquid  in  a  triangular 
groove  using  asymptotic  methods  and  a  regression  analysis.  Kolodziej  et  al.  (1999) 
used  the  boundary  collocation  method  to  solve  the  same  problem,  except  that  the 
liquid-vapor  interface  was  not  assumed  to  be  circular.  The  present  solution  was 
compared  to  that  obtained  by  Ayyaswamy  et  al.  for  9  =  5°  and  60°  for  the  full  range 
of  meniscus  contact  angle  (0.1°  <  <p  <  90°  -  9).  The  maximum  percent  differences 
for  9  =  5°  and  60°  were  3.7%  and  0.9%,  respectively.  In  comparison  to  the  results 
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by  Romero  and  Yost  (1996),  the  maximum  percent  differences  were  2.2%  for  0  =  5° 
and  2.3%  for  9  =  60°.  The  maximum  differences  in  the  comparison  with  the  results 
provided  by  Kolodziej  et  al.  for  the  lowest  value  of  Bond  number  presented  (Bo 
=  0.001)  was  4.5%  for  0  =  5°,  and  19.9%  for  9  =  60°.  This  may  be  due  to  the 
approximate  nature  of  the  solution  by  Kolodziej  et  al.,  which  was  in  terms  of  a 
truncated  infinite  summation. 

DiCola  (1968)  presented  the  solution  for  the  Poiseuille  number  for  the  laminar 
flow  of  a  constant  property  fluid  within  a  rectangular  groove.  While  interfacial  shear 
stress  at  the  liquid-vapor  interface  was  accounted  for,  the  groove  was  assumed  to  be 
completely  full,  with  a  meniscus  contact  angle  of  (f>  =  90°.  The  comparison  between 
the  equation  by  DiCola  and  the  results  of  the  present  analysis  for  laminar  flow  in  a 
family  of  rectangular  grooves  at  the  full  groove  condition  (</>  =  90°,  0.1  <  /3  <  1.0, 
T*  =  -0.1,  0.0,  and  1.0)  resulted  in  a  maximum  percent  difference  of  2.3%. 

The  present  model  was  compared  to  the  correlation  presented  by  Lin  and  Faghri 
(1997),  where  the  friction  factor  for  the  flow  of  liquid  in  triangular  grooves  with  liquid- 
vapor  shear  was  presented.  Unfortunately,  not  enough  information  was  provided  by 
Lin  and  Faghri  to  precisely  determine  the  limits  of  applicability  for  their  equation. 
Therefore,  the  correlation  was  evaluated  over  a  fairly  wide  range  for  comparison  with 
the  present  solution.  In  the  range  of  0.075  <  -t,*„  <  0.1,  the  maximum  percent 
difference  was  2.2%. 

1.4  Results  and  Discussion 
1.4.1  Parametric  Analysis 

A  numerical  study  has  been  completed  in  which  the  flow  field  in  a  trapezoidal 
groove  has  been  solved.  Specifically,  values  of  the  mean  velocity,  Poiseuille  number, 
and  volumetric  flow  rate  are  reported  for  various  values  of  the  groove  aspect  ratio, 
groove  half-angle,  meniscus  contact  angle,  and  dimensionless  shear  stress  at  the  liquid- 
vapor  interface  (Tables  8-25).  Figures  2  and  3  show  contour  plots  of  the  dimensionless 
velocity  obtained  for  two  trapezoidal  groove  geometries.  Figure  2  shows  that  the  flow 
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behavior  changes  significantly  with  shear  stress,  even  though  the  variation  of  rfv  is 
relatively  small  for  this  case.  The  maximum  velocity  is  located  along  the  liquid- 
vapor  interface  for  r*v  =  0  and  0.25,  and  within  the  interior  of  the  flow  field  for 
t*v  =  -0.25.  For  countercurrent  flow  [Fig.  2(c)],  a  region  of  reversed  flow  occurs  near 
the  intersection  of  the  groove  wall  and  the  liquid- vapor  interface.  Figure  3  shows 
the  velocity  field  in  which  the  fluid  depth  is  less  than  that  in  Fig.  2.  In  addition, 
the  shear  stress  t*v  varies  over  a  wider  range  in  Fig.  3.  In  contrast  to  Fig.  2,  the 
maximum  velocity  does  not  occur  at  the  centerline  of  the  groove,  but  near  the  corner 
of  the  groove  due  to  the  greater  depth  at  that  location.  The  length  of  the  groove  wall 
that  is  affected  by  the  reversed  flow  is  greater  in  Fig.  3(c)  due  to  the  fact  that  the 
angle  between  the  liquid-vapor  interface  and  the  groove  wall  is  more  acute  than  that 
in  Fig.  2(c). 

Figure  4  and  Tables  8,  11,  and  14  present  the  mean  velocity  versus  shear  stress  at 
the  liquid- vapor  interface  for  several  values  of  the  groove  half-angle.  The  range  of  the 
meniscus  contact  angle  (0  <  (j)  <  90°  -  6)  was  divided  equally  to  show  the  behavior 
of  the  mean  velocity  with  4 t>.  The  mean  velocity  increases  linearly  with  shear  stress 
since  the  flow  is  aided  by  t*v.  In  addition,  v*  increases  with  groove  half-angle  and 
meniscus  contact  angle,  which  is  a  result  of  an  increase  in  cross-sectional  area.  As 
the  groove  half-angle  6  increases,  v*  becomes  more  sensitive  to  the  meniscus  contact 
angle  <j>  due  to  the  increased  length  of  the  perimeter  of  the  liquid-vapor  interface. 

Figure  5  and  Tables  9,  12,  and  15  show  that  the  Poiseuille  number  decreases 
dramatically  as  the  shear  stress  at  the  liquid-vapor  interface  increases.  For  counter- 
current  flow  (t*v  <  0),  the  Poiseuille  number  is  a  strong  function  of  shear  stress  since 
the  mean  velocity  approaches  zero.  In  addition,  the  Poiseuille  number  decreases  with 
increasing  meniscus  contact  angle  for  a  given  value  of  shear  stress.  For  cocurrent  flow 
(r,*  >  0),  the  Poiseuille  number  is  a  lesser  function  of  the  shear  stress,  but  increases 
significantly  with  meniscus  contact  angle.  The  Poiseuille  number  is  a  weak  function 
of  the  groove  half-angle. 

The  volumetric  flow  rate  versus  shear  stress  for  various  meniscus  contact  angles 
can  be  seen  in  Fig.  6  and  Tables  10,  13,  and  16.  The  volumetric  flow  rate  and  mean 
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Figure  2:  Dimensionless  velocity  fields  for  laminar  flow  in  trapezoidal  grooves  (/? 
1.0,  (j)  =  40°,  6  =  30° ):  (a)  r,*  =  0.25  (cocurrent  flow);  (5)  t*v  =  0.0;  (c)  Tj*  =  -0 
(countercurrent  flow). 


Figure  3:  Dimensionless  velocity  fields 
1.0,  d>  =  10°,  6  =  45°);  (a)  r*  =  5.0  ( 
(countercurrent  flow). 


laminar  flow  in  trapezoidal  grooves  (/? 
irrent  flow);  (6)  t*v  =  0.0;  (c)  t*v  -  -i 


velocity  display  similar  trends.  The  volumetric  flow  rate  is  a  linear  function  of  shear 
stress,  and  increases  with  meniscus  contact  angle  and  groove  half-angle.  The  flow 
rate  is  slightly  more  sensitive  than  the  mean  velocity  with  respect  to  the  meniscus 
contact  angle. 

The  mean  velocity,  Poiseuille  number  and  volumetric  flow  rate  are  given  as  func¬ 
tions  of  the  interfacial  shear  stress  for  a  constant  meniscus  contact  angle  in  Figs.  7, 
8,  and  9  and  in  Tables  17-25.  For  0  =  0  and  30°,  the  mean  velocity  increases  and 
then  decreases  with  groove  aspect  ratio.  This  point  will  be  further  elucidated  in  fol¬ 
lowing  figures.  The  Poiseuille  number  (Fig.  8)  is  a  weak  function  of  both  0  and  (3  for 
cocurrent  flow,  but  varies  significantly  with  (3  for  countercurrent  flow  in  rectangular 
grooves  (0  =  0°).  The  volumetric  flow  rate  (Fig.  9)  follows  the  same  trends  as  the 
mean  velocity,  but  is  more  heavily  influenced  by  /3,  and  0. 

Figure  10  and  Table  28  present  the  mean  velocity,  Poiseuille  number  and  volumet¬ 
ric  flow  rate  for  cj>  =  30°  and  r*v  =  5.0.  As  mentioned  previously,  the  mean  velocity 
increases  and  then  decreases  with  (3  for  0  <  30°.  This  phenomenon  also  impacts  the 
volumetric  flow  rate,  where  V*  attains  a  maximum  value  with  respect  to  6  for  a  given 
groove  half-angle  <  15°.  Figure  11  and  Table  29  show  v*,  Po,  and  V*  versus  0  for 
(3  =  1.0  and  t*v  =  5.0.  In  general,  these  functions  increase  with  both  meniscus  contact 
angle  and  groove  half-angle,  except  for  the  Poiseuille  number  for  0  =  0°  and  cj)  <20°. 

1.4.2  Semi-Analytical  and  Two-Point  Numerical  Solutions  for  v* 

As  seen  in  Figs.  4  and  7,  the  mean  velocity  is  a  linear  function  of  the  imposed  shear 
stress  at  the  liquid-vapor  interface.  Since  a  direct  numerical  simulation  of  the  liquid 
flow  field  for  a  number  of  values  of  the  shear  stress  is  computer  resource  intensive, 
it  is  appropriate  to  seek  a  semi-analytical  expression  for  v*.  Figure  12(a)  shows  the 
definition  of  the  parameters  involved,  where  the  mean  velocity  when  the  shear  stress 
is  zero  (i^)  is  given  by  the  numerical  solution.  The  value  for  the  liquid- vapor  shear  for 
which  the  mean  velocity  is  zero  (r*v  0)  is  given  by  the  following  force  balance  analysis. 
Figure  12(6)  shows  a  differential  element  of  the  liquid  in  the  groove.  A  force  balance 
between  the  pressure  drop  and  the  shear  forces  at  the  liquid-vapor  interface  and  at 
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Figure  5:  Po  versus  rt*v  for  laminar  flow  in  trapezoidal  grooves  (/?  =  1.0):  (o)  8  —  O' 
(i b )  8  =  30°;  (c)  8  =  60°. 


Figure  9:  V 
(1 b )  6  =  30°; 


-1  0 


1 


' lv 


3 


4  5 


*  versus  r;*  for  laminar  flow  in  trapezoidal  grooves  (<p  =  30° ):  (a)  9  =  0° 
(c)  6  =  60°. 
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the  wall  results  in  the  following  relation. 


PyA[  Py+dyAl  +  TlvAlv  T^jAyj  0  (14) 

The  areas  over  which  the  shear  stresses  rlv  and  tT  act  are  Alv  =  Ptvdy  and  Aw  =  Pdy, 
respectively.  Using  these  areas  and  nondimensionalizing  gives 

a;  +  -  kp'  =  o  (is) 

For  Poiseuille  flow  in  ducts  of  arbitrary  cross  section,  and  combined  Couette-Poiseuille 
flow  between  flat  plates,  the  shear  stress  at  the  wall  is  related  to  the  mean  velocity 
of  the  fluid  by  a  constant  (White,  1991).  In  the  present  analysis,  it  is  assumed  that 
this  also  holds  for  the  flow  of  liquid  in  a  trapezoidal  groove  with  an  imposed  shear 
stress  at  the  liquid- vapor  interface. 

=  Cyv*  (16) 


It  should  be  noted  that  the  constant  Ci  is  probably  a  function  of  the  groove  geom¬ 
etry  and  meniscus  contact  angle.  However,  since  the  objective  of  this  analysis  is  to 
determine  the  liquid-vapor  shear  stress  when  the  mean  liquid  velocity  is  zero,  this 
functionality  is  unimportant. 

Substituting  this  relation  into  the  force  balance  equation  results  in  the  following 
expression  for  mean  velocity 

F  =  cTF  (A‘  +  T,;F‘;)  <17) 


(18) 


where  the  perimeter  of  the  liquid- vapor  interface  is 

2i?*sin-1  )  for  6  +  <j>  <  tt/2 

2  (/3  +  tan0)  for  6  +  <j>  =  7r/2 
The  mean  velocity  is  zero  when  the  shear  stress  at  the  liquid- vapor  interface  is 

A ‘  (19) 


Tlv,  0  — 


tv 
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Figures  13(a)  and  13(6)  show  the  results  of  eqn.  (19).  The  numerical  results  shown 
in  Figs.  4  and  7  were  extrapolated  to  determine  the  values  for  shear  stress  at  the 
liquid-vapor  interface  when  v*  =  0.  The  prediction  given  by  eqn.  (19)  is  quite  good 
given  the  simplicity  of  the  closed-form  solution.  The  equation  for  the  normalized 
mean  velocity  as  a  function  of  the  shear  stress  is  given  by 

v'  =  v*/vq  =  1  —  t'  (20) 


where  r'  =  t*v/t*v0.  The  semi-analytical  solution  for  the  normalized  mean  velocity 
is  shown  in  Fig.  13(c)  with  the  corresponding  numerical  data  presented  in  Figs.  4 
and  7.  Equation  (20)  predicts  93%  of  the  data  to  within  ±30%  over  the  range  of  the 
meniscus  contact  angle,  groove  half-angle,  groove  aspect  ratio  and  liquid-vapor  shear 
stress  examined  in  Figs.  4  and  7. 

The  two-point  numerical  solution  of  v*  as  a  function  of  t*v  is  also  shown  in  Fig. 
12(a).  This  involves  computing  the  flow  field  for  two  values  of  liquid- vapor  shear 
stress.  The  equation  for  the  normalized  mean  velocity  using  this  solution  is  given  by 
eqn.  (20),  but  the  shear  stress  when  the  mean  velocity  is  zero  is 


lv,°  (i  -  KK) 


(21) 


1.4.3  Regression  Analysis 

A  linear  regression  analysis  was  performed  to  determine  the  mean  velocity  as 
a  function  of  the  groove  geometry,  meniscus  contact  angle,  and  shear  stress  at  the 
liquid-vapor  interface.  Equation  (22)  and  Table  1  show  the  results  of  this  analysis. 
The  regression  was  found  by  determining  the  relationships  between  the  mean  velocity 
and  each  of  the  groove  aspect  ratio,  groove  half  angle,  meniscus  contact  angle,  and 
interfacial  shear  stress  separately,  then  summing  these  relations  to  arrive  at  the  model 
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Figure  13:  Numerical  and  semi-analytical  solution  comparison:  (a)  -r,*  >0  versus  0  for 
=  1.0  (from  Fig.  4);  (6)  -7^>0  versus  /3  for  0  =  30°  (from  Fig.  7);  (c)  Normalized 
mean  velocity  versus  normalized  shear  stress  at  the  liquid- vapor  interface  (from  Figs. 
4  and  7). 
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Figure  14:  Comparison  of  data  with  regression, 
shown  in  eqn.  (22). 

v*  =  adz  +  b92  +  c92(3  +  d92<f>  +  e92r*v  +  fd<f>  +  gO  +  h(f>3  +  i<p2d  +  + 

k<f>rfv  +  l<p9r*v  +  m<t>  +  n(32  +  op29  +  +  qpO  +  rPr{v  +  (22) 

s /36(f)  +  tpdT*v  +  up(f)T*v  +v(3  +  WTIV  +  xp0(pT*v  +  y 

A  statistical  software  package  SPSS  (Statistical  Package  for  the  Social  Sciences)  was 
utilized  to  perform  the  regression  analysis,  and  the  outcome  for  the  coefficients  is 
shown  in  Table  1.  A  comparison  of  the  observed  mean  velocity  versus  the  predicted 
mean  velocity  is  presented  in  Fig.  14.  The  regression  predicts  85.0%  of  the  observed 
data  presented  in  Figs.  4,  7,  10(a),  and  11(a)  to  within  ±30%. 

1.4.4  Effect  of  Groove  Fill  Ratio 

Figure  15(a)  shows  the  case  when  liquid  evaporates  from  a  trapezoidal  groove. 
Initially,  the  groove  is  full  with  (f>  +  9  =  90°.  The  contact  angle  decreases  until 
the  minimum  meniscus  contact  angle  </>o  for  the  particular  solid-liquid  combination  is 
reached.  Past  this  point,  the  meniscus  detaches  from  the  top  of  the  groove  and  recedes 
until  the  lowest  part  of  the  meniscus  reaches  the  bottom  of  the  groove  (Hopkins  et 
al.,  1999).  When  the  thickness  of  the  liquid  film  at  the  bottom  of  the  groove  is 
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Table  1:  Coefficients  for  regression  analysis. 


Coefficient 

Value 

Coefficient 

Value 

a 

4.045E-6 

m 

2.492E-3 

b 

-4.747E-4 

n 

-0.1780 

c 

5.656E-5 

0 

3.186E-3 

d 

3.540E-6 

p 

-9.241E-2 

e 

3.970E-5 

Q 

-1.031E-2 

f 

-2.299E-4 

r 

7.430E-2 

9 

K WEM 

s 

8.072E-5 

h 

-9.312E-8 

t 

-9.775E-4 

i 

7.011E-7 

u 

4.309E-3 

j 

-1.495E-5 

V 

0.3420 

k 

-1.460E-3 

w 

8.266E-2 

i 

3.194E-5 

X 

4.224E-5 

y 

-0.1780 

on  the  order  of  several  hundred  Angstroms,  forces  due  to  London- van  der  Waals 
interactions  with  the  surrounding  liquid  and  solid  molecules  induce  instabilities  in 
the  fluid  (Ruckenstein  and  Jain,  1974).  These  instabilities  cause  the  liquid  in  the 
groove  to  bifurcate  into  two  separate  flows  in  the  corners  of  the  groove,  which  are 
each  equivalent  to  the  flow  in  a  triangular  groove.  The  liquid  in  the  two  corners  of 
the  groove  will  continue  to  recede  until  it  is  depleted. 

The  dimensions  of  the  grooves  analyzed  by  Castle  et  al.  (2000)  were  used  to 
determine  the  volumetric  flow  rate  of  ethanol  in  a  trapezoidal  copper  groove  as  a 
function  of  the  amount  of  liquid  in  the  groove.  Faghri  (1995)  gives  ©0  =  7°  for  a 
receding  meniscus  of  ethanol  on  copper.  In  terms  of  the  present  analysis,  as  the 
liquid  recedes  into  the  groove,  the  groove  aspect  ratio  j5  increases.  In  addition,  after 
bifurcation  occurs  /3  =  0,  and  the  groove  half-angle  6  changes  to  that  of  the  corner  of 
the  groove.  A  relation  for  the  point  at  which  bifurcation  occurs  is  provided  where  the 
liquid  is  assumed  to  bifurcate  when  the  lowest  part  of  the  meniscus  actually  reaches 
the  bottom  of  the  groove.  The  radius  of  curvature  at  the  bifurcation  point  for  the 
trapezoidal  groove  is 

R*b  =  -  [l  +  (/?  +  tan  0)2]  (23) 
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Figure  15:  Effect  of  groove  fill  ratio  on  liquid  flowing  in  a  trapezoidal  groove:  (a) 
Parametric  analysis  (to  scale,  h  =  0.03831  cm,  w  =  0.03445  cm,  6  —  14.62°);  ( b ) 
Definition  of  variables  prior  to  bifurcation  of  the  liquid;  (c)  Definition  of  variables 
after  bifurcation  of  the  liquid. 


28 


Figure  16:  Meniscus  contact  angle  (in  arc  degrees)  at  bifurcation  versus  #  for  various 
values  of  (3. 

The  meniscus  contact  angle  at  bifurcation  as  a  function  of  the  groove  geometry  is 

-l 


06  =  tan 


-l 


sin# 


cos#  + 


2  sin  9((3  +  tan  #) 

1  —  (/5  +  tan  #)2  J 


-l 


cos# 


(24) 


=  tan 


-l 


(3 (32  —  4)  cos  #  +  01  cos  3# 


|_8/?cos#+  (/?2  +  4)  sin#  +  /?2  sin  3#  J 
This  function  is  shown  in  Fig.  16  for  various  values  of  (3.  For  (3  <  1  (slender 
grooves),  the  contact  angle  is  negative  for  all  values  of  #,  which  indicates  that  the 
lowest  part  of  the  meniscus  cannot  reach  the  groove  bottom  and  still  maintain  a 
circular  shape.  For  (3  =  1,  the  only  value  of  0&  for  which  the  meniscus  will  reach  the 
groove  bottom  is  #  =  0.  If  the  groove  aspect  ratio  is  /?  >  1,  the  meniscus  can  reach 
the  groove  bottom  over  a  range  of  #  values. 

The  dimensional  height  and  width  of  the  liquid  before  and  after  bifurcation  are 
defined  in  Figs.  15(6)  and  15(c).  For  a  given  liquid  height,  the  width  of  the  liquid  is 


Wi 


(25) 


w  +  2 hi  tan  #  before  bifurcation 
2 hi  tan  #  after  bifurcation 

The  radius  of  curvature  of  the  liquid- vapor  interface  and  the  cross-sectional  area  of 
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the  liquid  are 


R=< 


Wi 

T 


1  +  <  cot  #  — 


sin0 


[sin  #sin(#  +  <fi) 


]} 


-2 


before  bifurcation 


(26) 


hi  tan0W  1  +  <  cot  6  — 


{■ 


sin  (j) 


[sin  #  sin  (#  +  4>) 


]}' 


after  bifurcation 


|  (10  +  w,)  -  R2  cos-1  ^1  -  (H)  +  ^  \/l  -  (^j 


2R; 

before  bifurcation 


hf  tan  6  —  R2  cos  1  \  1  — 


'hi  tan#' 

J?  . 


+  hiRtan9\  l  — 


'hi  tan#' 

i? 


after  bifurcation 


(27) 


The  cross-sectional  area  of  the  groove  is 


Ag  =  h(w  +  h  tan  #) 


(28) 


Table  30  gives  the  geometric  values  of  the  parametric  analysis  depicted  in  Fig.  15(a). 

The  radius  of  curvature  of  the  liquid- vapor  interface  is  shown  in  Fig.  17(a)  as 
a  function  of  the  groove  fill  ratio.  When  the  groove  is  nearly  full,  the  radius  of 
curvature  approaches  infinity.  As  the  amount  of  liquid  in  the  groove  decreases,  the 
radius  of  curvature  is  relatively  constant,  and  then  becomes  very  small  after  the  liquid 
bifurcates  into  the  corners  of  the  groove.  Figures  17(5)  and  17(c)  and  Table  31  show 
two  flow  parameters,  which  were  defined  as  follows. 


^1  = 


fJLV 

(- dp/dy ) 


(29) 


(-dp/dy) 


(30) 


These  parameters  allow  the  presentation  of  the  mean  velocity  and  volumetric  flow 
rate  before  and  after  bifurcation  on  the  same  graph.  Both  flow  parameters  increase 
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monotonically  with  groove  fill  ratio,  as  expected.  In  Fig.  17(c),  for  a  groove  fill  ratio 
of  AiJAg  =  0.158,  the  volumetric  flow  rate  was  1%  of  that  for  the  full  groove  due  to 
the  decrease  in  flow  area.  This  figure  shows  that  the  groove  was  essentially  shut  off  for 
Ai/Ag  <  0.158,  which  was  just  prior  to  bifurcation.  Using  the  functional  relationship 
for  mean  velocity  given  by  eqn.  (17),  the  volumetric  flow  rate  is 

V*  =  AL  (a;  +  t;vp;v)  pi) 

When  t*v  =  0,  V  oc  A],  which  confirms  that  the  volumetric  flow  rate  should  decrease 
rapidly  with  decreasing  flow  area  for  Ai  <  1. 


1.4.5  Capillary  Limit  Analysis  for  a  Revolving  Helically- Grooved  Heat 
Pipe 

Using  the  results  of  the  numerical  analysis,  the  capillary  limit  prediction  for  a  re¬ 
volving  helically-grooved  heat  pipe  proposed  by  Thomas  et  al.  (1998)  was  improved 
by  accounting  for  the  effects  of  working  fluid  fill  amount  and  the  shear  stress  at  the 
liquid-vapor  interface.  The  improved  model  was  compared  to  the  experimental  data 
collected  by  Castle  et  al.  (2000),  who  determined  the  capillary  limit  of  a  revolving 
helically-grooved  copper-ethanol  heat  pipe  for  radial  accelerations  of  |c£|  =  0.01,  2.0, 
4.0,  6.0,  8.0  and  10.0-g  and  groove  fill  ratios  of  Vi/Vg  =  0.5,  1.0  and  1.5.  The  liquid 
inventory  volume  and  the  total  groove  volume  are  given  by  Vi  and  Vg,  respectively. 
The  dimensions  of  the  heat  pipe  examined  by  Castle  et  al.  (2000)  are  given  in  Table 
2.  A  pressure  balance  within  the  heat  pipe  results  in  the  following  expression  for  the 
capillary  limit  (Faghri,  1995;  Chi,  1976) 


Apcap.max  ^  Apr  H"  Ap;  "k  Apbf 


(32) 


The  maximum  capillary  pressure  for  an  axial  groove  is 


ApCap,max  —  _ 
JXr 


(33) 


It  is  assumed  that  the  capillary  limit  occurs  when  the  liquid  bifurcates  into  the  corners 
of  the  grooves.  This  statement  is  based  on  the  results  of  the  variation  of  the  volumetric 
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Figure  17:  Effect  of  groove  fill  ratio:  (a)  Radius  of  curvature  of  the  liquid-vapor 
interface;  (6)  Mean  velocity  parameter;  (c)  Volumetric  flow  rate  parameter. 


Table  2:  Specifications  of  the  heat  pipe  test  article  examined  by  Castle  et  al.  (2000). 


Groove  Height,  h 

0.03831  cm 

Groove  Base  Width,  w 

0.03445  cm 

Groove  Half-Angle,  9 

14.62° 

Evaporator  Length,  Le 

15.2  cm 

Adiabatic  Length,  La 

8.2  cm 

Condenser  Length,  Lc 

15.2  cm 

Vapor  Space  Radius,  Rv 

0.6795  cm 

Number  of  Grooves,  Ng 

50 

Helical  Pitch  Length,  Lp 

135.8  cm 

Helix  Radius,  Rh 

0.6992  cm 

flow  rate  with  groove  fill  ratio  in  Fig.  17.  The  capillary  radius  when  the  liquid 
bifurcates  is  given  by  eqn.  (23). 

Rc  =  —  [l  +  (P  +  tan  0)2]  (34) 


The  meniscus  contact  angle  at  bifurcation  is  a  function  of  the  groove  geometry,  as 
presented  in  eqn.  (24).  These  two  equations  were  solved  to  determine  the  capillary 
radius  at  the  point  of  dryout  in  the  following  manner.  The  groove  aspect  ratio 
was  calculated  using  the  groove  height  and  width  found  in  Table  2  (/?  =  0.4496). 
This  value  of  f3  indicates  that  the  meniscus  must  recede  into  the  groove  before  the 
bifurcation  point  can  be  reached  (Fig.  16).  Therefore,  eqn.  (24)  was  solved  iteratively 
to  determine  the  aspect  ratio  for  which  the  meniscus  contact  angle  was  equal  to  the 
minimum  meniscus  contact  angle  for  the  copper-ethanol  system.  The  resulting  aspect 
ratio  was  (3  —  1.2112.  Solving  for  the  height  of  the  meniscus  in  the  groove  results  in 
h  =  hi  =  0.01422  cm.  Substituting  this  value  into  eqn.  (34)  for  the  meniscus  radius 


of  curvature  gives  Rc  =  0.02252  cm. 

For  a  circular  cross  section  heat  pipe  assuming  laminar,  incompressible  flow  with 
uniform  heat  input  and  output  along  the  lengths  of  the  evaporator  and  condenser, 
respectively,  the  pressure  drop  in  the  vapor  is 


SflyLesQt 

7T  pyhfgRy 


(35) 
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The  pressure  drop  in  the  liquid  was  found  by  using  the  normalized  mean  velocity 
relation  [eqn.  (20)]  rewritten  in  dimensional  form. 


7  |  PlVl  ,  ^IV  \  Ant 

dpi  =  ~  \  Tr=z  +  TT7-  I  dV 


(36) 


hfv 5  htTfv  |0 

In  a  heat  pipe,  the  liquid  flows  opposite  to  the  vapor  in  all  regions.  Therefore,  the 
shear  stress  at  the  liquid- vapor  interface  for  countercurrent  flow  was  used. 

Pv  fa? 


7~lv  — 


fv 


(37) 


The  Poiseuille  number  of  the  vapor  flow  was  modeled  as  laminar  flow  within  a  smooth 
tube  with  a  circular  cross  section 


Po„  =  fvRev  =  16 


(38) 


Substituting  these  relations  into  eqn.  (36)  gives 

PlVl  4 


dpi  =  - 


dy 


(39) 


hfv  q 

It  is  assumed  that  the  cross-sectional  area  of  the  liquid  is  constant  along  the  length  of 
the  groove.  For  a  constant  heat  flux  in  both  the  evaporator  and  condenser  sections, 
eqn.  (39)  can  be  integrated  to  determine  the  total  pressure  drop  in  the  helical  groove. 


a  t  I  W^i»niax  4/i  v^t;,max 

A Pl  =  ^eff  72= - I  p  -* 

hfvQ  hiKvTiv0 


m  *> 


(40) 


The  maximum  liquid  velocity  in  a  groove  is 


^,max  — 


Similarly,  the  maximum  vapor  velocity  is 


Qg 


PlAlhfe 


Qt 


(41) 


(42) 


t-'t?,max  7*>o  l 

7lR^pvhfg 

where  the  total  heat  transported  by  the  heat  pipe  accounts  for  the  contributions  by 
all  of  the  individual  grooves. 
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Using  the  above  relations,  the  liquid  pressure  drop  in  a  groove  as  a  function  of  the 
transported  heat  is 


Ap;  = 


Le{(  (  hlQg _ 4  jdyQt 


‘lv,  0 


(44) 


hfg  yhfv*0piAi  nhiRlpv'i 

The  body  forces  imposed  on  the  fluid  within  a  particular  groove  may  either  aid 
or  hinder  the  return  of  the  fluid  to  the  evaporator,  depending  on  the  groove  pitch  Lp 
and  the  circumferential  location  of  the  starting  point  of  the  helical  groove  (Thomas  et 
ah,  1998;  Klasing  et  ah,  1999).  However,  even  if  the  body  force  hinders  the  return  of 
the  fluid,  each  groove  contributes  to  the  heat  transported  Qt ■  Therefore,  the  capillary 
limit  equation  [eqn.  (32)]  was  first  solved  for  the  heat  transported  by  each  individual 
groove  Qg ,  and  the  results  were  summed  to  determine  the  total  heat  transport  Qt- 
Since  the  pressure  drop  in  the  vapor  space  and  the  pressure  drop  in  each  groove  were 
functions  of  the  total  heat  transport,  eqn.  (32)  was  solved  iteratively.  The  body 
forces  due  to  acceleration  and  gravity  were  integrated  over  the  length  of  the  groove 
to  find  the  average  pressure  drop  (Thomas  et  al.,  1998). 
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ds 
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Combining  the  above  relations,  the  general  expression  for  the  maximum  capillary 
limit  for  a  single  helical  groove  which  accounts  for  shear  stress  at  the  liquid-vapor 
interface  and  the  effect  of  groove  fill  ratio  is  given  by 
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A  closed-form  solution  for  the  capillary  limit  of  a  heat  pipe  with  straight  axial 
grooves  and  no  body  forces  can  be  derived  from  eqn.  (46). 
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Figure  18(a)  shows  the  results  of  the  closed-form  solution  [eqn.  (47)].  The  groove 
geometry  given  by  Castle  et  al.  (2000)  was  used,  except  that  straight  axial  grooves 
were  assumed  ( Lv  — >  oo)  instead  of  helical  grooves.  Table  32  gives  the  results  of 
the  computer  model  for  the  geometry  analyzed  by  Castle  et  al.  (2000)  for  a  shear 
stress  of  r*  =  5.0.  This  information  was  used  as  part  of  the  two-point  numerical 
solution,  along  with  the  data  presented  in  Table  31  for  t>J.  Over  the  range  of  groove 
fill  ratio  examined,  the  capillary  limit  increased  with  Ai/Ag  by  more  than  three  orders 
of  magnitude.  For  this  case,  the  semi-analytical  solution  and  the  two-point  numerical 
solution  were  nearly  identical  due  to  the  low  vapor  velocities,  and  hence  the  low 
liquid-vapor  shear  stress.  This  point  is  further  demonstrated  in  Fig.  18(a)  by  the 
graph  indicated  by  “No  Shear,”  where  the  term  in  eqn.  (47)  that  accounts  for  the 
effect  of  liquid-vapor  shear  on  the  liquid  was  dropped  by  allowing  t*v  0  — >•  — oo.  The 
maximum  difference  in  the  closed  form  solution  compared  to  the  two-point  numerical 
solution  was  5%.  The  capillary  limit  decreases  when  shear  stress  is  accounted  for,  as 
expected.  Figure  18(6)  shows  the  closed-form  solution  when  water  is  the  working  fluid. 
In  this  case,  the  difference  between  the  results  of  the  two-point  numerical  solution 
and  the  no-shear  solution  is  much  more  pronounced  (45%  maximum  difference)  due 
to  the  significantly  higher  vapor  velocities  involved.  The  agreement  between  the 
semi-analytical  solution  and  the  two-point  numerical  solution  is  quite  good.  The 
semi-analytical  solution  offers  very  close  results  with  a  significantly  reduced  amount 
of  computer  resources  required. 

The  capillary  limit  prediction  for  the  helically-grooved  heat  pipe  given  by  eqn. 
(46)  is  shown  in  Fig.  19  using  the  semi-analytical  solution  [eqn.  (19)]  for  |ar|  =  0.0 
and  10.0-g.  The  capillary  limit  heat  transfer  increases  both  with  groove  fill  ratio 
and  working  temperature.  In  addition,  Qcap  increases  significantly  with  the  radial 
acceleration  due  to  the  improved  liquid  pumping  ability  of  the  helical  grooves  (Thomas 
et  al.,  1998). 

The  experimental  data  collected  by  Castle  et  al.  (2000)  for  the  capillary  limit 
of  a  revolving  helically-grooved  heat  pipe  versus  radial  acceleration  rates  are  shown 
in  Fig.  20,  along  with  the  predictions  of  the  present  semi-analytical  model  and  that 
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Figure  18:  Maximum  heat  transport  predicted  by  the  closed-form  solution  versus 
groove  fill  ratio  (Straight  axial  grooves,  no  body  forces,  Tsat  =  40°C):  (a)  Ethanol; 
(b)  Water. 
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Figure  19:  Maximum  heat  transport  versus  groove  fill  ratio  for  several  working  tem 
peratures  (Ethanol):  (a)  |ar|  =  0.0-g;  ( b )  |or|  =  10.0-g. 


given  by  Castle  et  al.  (2000).  During  the  experiments,  the  working  temperature  was 
not  held  constant,  so  the  present  model  was  evaluated  at  the  saturation  temperature 
reported  for  a  given  value  of  radial  acceleration.  For  a  groove  fill  ratio  of  AjAg  =  0.5 
[Fig.  20(a)],  the  present  model  more  closely  matches  the  experimental  data  than 
the  model  by  Castle  et  al.  (2000),  which  did  not  account  for  the  groove  fill  ratio  or 
liquid- vapor  shear  stress.  The  maximum  difference  in  the  experimental  and  analytical 
results  by  Castle  et  al.  (2000)  was  418%,  whereas  the  maximum  difference  for  the 
current  model  was  104%.  For  A[/Ag  =  1.0  [Fig.  20(5)],  the  present  model  overpredicts 
the  experimental  data  by  a  maximum  difference  of  158%  and  the  previous  model  by 
87%. 

1.5  Conclusions 

A  numerical  study  has  been  concluded  where  the  mean  velocity,  Poiseuille  number, 
and  volumetric  flow  rate  of  liquid  in  a  trapezoidal  groove  have  been  determined  as 
functions  of  groove  geometry,  meniscus  contact  angle  and  shear  stress  at  the  liquid- 
vapor  interface.  The  mean  velocity  and  volumetric  flow  rate  have  been  shown  to 
be  linear  functions  of  shear  stress,  and  the  Poiseuille  number  is  a  strong  function 
of  the  shear  stress  for  countercurrent  flow.  A  semi-analytical  solution  and  a  two- 
point  numerical  solution  for  the  mean  velocity  were  presented  and  used  to  predict 
the  capillary  limit  of  a  revolving  helically-grooved  heat  pipe  for  various  groove  fill 
ratios.  Interfacial  shear  stress  due  to  countercurrent  flow  in  a  heat  pipe  decreases  the 
maximum  heat  transport.  For  cases  in  which  the  vapor  velocities  are  high,  this  effect 
is  more  pronounced.  The  groove  fill  ratio  was  shown  to  have  a  significant  impact  on 
heat  pipe  performance.  Underfilling  the  heat  pipe  examined  by  10%  resulted  in  a 
decrease  in  the  predicted  capillary  limit  by  approximately  17  to  20%  for  water  and 
ethanol,  respectively. 
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Figure  20:  Maximum  heat  transport  versus  radial  acceleration:  (a)  At/Ag  =  0.5;  (6) 
At/Ag  =  1.0. 
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2  Fully  Developed  Laminar  Flow  in  Sinusoidal  Grooves 
2.1  Introduction 

Isotropic  materials  such  as  quartz  glass  or  borosilicate  glass  can  be  chemically 
etched  to  form  micro-grooves  for  the  enhancement  of  evaporative  heat  transfer  in 
chip-level  devices  (Kirshberg  et  al.,  1999).  A  typical  etched  profile  in  glass  is  shown 
in  Fig.  21,  where  the  profile  is  smooth  instead  of  having  sharp  corners  seen  in  the 
etching  of  crystalline  materials  such  as  silicon  (Maluf,  2000).  While  many  studies  have 
been  performed  on  sharp-cornered  geometries  [triangular  grooves  (Ayyaswamy  et  al., 
1974;  Ma  et  al.,  1994;  Romero  and  Yost,  1996;  Lin  and  Faghri,  1997;  Kolodziej  et  al., 
1999),  rectangular  grooves  (DiCola,  1968;  Schneider  and  DeVos,  1980;  Khrustalev 
and  Faghri,  1999),  and  trapezoidal  grooves  (Hopkins  et  al.,  1999)],  very  little  in¬ 
formation  is  available  in  the  open  literature  on  the  flow  of  liquid  in  rounded-corner 
geometries.  Stroes  and  Catton  (1997)  compared  the  capillary  performance  of  trian¬ 
gular  and  sinusoidal  grooves  by  means  of  an  experimental  study.  Two  sets  of  grooves 
were  machined  into  stainless  steel  test  plates  such  that  the  cross-sectional  areas  of 
the  grooves  were  equal.  Strip  heaters  were  placed  under  the  plates  to  provide  heat 
input.  The  test  plates  were  placed  at  inclination  angles  of  4  and  6°  and  ethanol  was 
added  to  the  grooves  until  the  liquid  reached  the  lands  of  the  grooves.  The  average 
wetted  length  of  each  set  of  grooves  was  recorded  as  the  heat  input  was  varied  from 
0  to  25  W.  The  study  showed  that  the  triangular  grooves  had  a  greater  capillary 
pumping  ability  compared  to  the  sinusoidal  grooves  with  the  same  cross-sectional 
area,  inclination  angle,  and  heat  input.  Stroes  and  Catton  postulated  that  this  was 
due  to  the  axial  rate  of  change  of  the  radius  of  curvature  of  the  meniscus.  Sinusoidal 
grooves,  however,  could  dissipate  a  given  heat  input  with  a  smaller  wetted  area  than 
triangular  grooves  due  to  the  larger  wetted  perimeter  found  with  sinusoidal  grooves. 
The  objective  of  the  present  study  was  to  examine  the  fully-developed  laminar  flow 
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Figure  21:  Grooves  chemically  etched  in  glass  (Courtesy  of  D.  Liepmann,  University 
of  California  at  Berkeley). 

of  liquid  in  sinusoidal  grooves.  The  effects  of  countercurrent  and  cocurrent  vapor 
flow  over  the  liquid-vapor  interface  were  investigated  by  relating  the  liquid  velocity 
gradient  to  the  friction  factor  of  the  vapor.  The  variation  of  the  shear  stress  on  the 
liquid-vapor  interface  (Khrustalev  and  Faghri,  1999)  was  neglected,  and  the  liquid- 
vapor  interface  was  assumed  to  be  circular  (Bo  1).  The  mean  velocity,  volumetric 
flow  rate  and  Poiseuille  number  were  determined  as  functions  of  the  interfacial  shear 
stress,  the  meniscus  contact  angle,  the  groove  aspect  ratio  and  the  amount  that  the 
groove  was  filled. 

2.2  Mathematical  Model 


A  constant  property  liquid  flows  steadily  in  a  sinusoidal  groove  as  shown  in  Fig. 
22.  A  meniscus,  which  is  assumed  to  be  circular,  comprises  the  liquid-vapor  interface. 
For  fully  developed  laminar  flow,  the  conservation  of  mass  and  momentum  equations 
reduce  to  the  classic  Poisson  equation  in  dimensionless  form  (White,  1991) 

d2v*  d2v* 


dx *2  dz*2 

On  the  groove  wall,  the  no-slip  condition  is  in  effect. 

v*  =  0  :  0<x*<w*/2,  z*  =  i  j  1  +  cos 


~n[J  +  1 


At  the  line  of  symmetry,  the  velocity  gradient  is  zero  in  the  x*  direction 
dv* 


dx * 


=  0:  x*  =  0,  0<z*  <h;  +R*\  1- 


2  R* 


-R* 


(48) 
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The  dimensionless  radius  of  curvature  is  given  by 

"■-(f) 

where 


(1  +  d*  tan  <f> )2 
(i d *  —  tan<£)2 


(51) 
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—  + 1 
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(52) 


The  maximum  value  for  the  meniscus  contact  angle  <j>  for  a  wetting  fluid  can  be 
determined  for  a  given  geometry  by  allowing  the  radius  of  curvature  to  approach 
R*  oo. 


0max  —  tan  ^  2/? 


,  wl 


(53) 


At  the  liquid-vapor  interface,  a  uniform  shear  stress  is  imposed  in  the  y  direction. 

=  r,; :  o  <  I*  <  to, 72, 2*  =  K  +  Jryi-Qfi)  -  ^  ~  **2  (54> 
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The  dimensional  shear  stress  at  the  liquid-vapor  interface  can  be  cast  in  terms  of  the 
friction  factor  of  the  vapor. 

Pv  (*T)2 
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fv  for  cocurrent  flow 


fv  for  countercurrent  flow 


(55) 


The  Poiseuille  number  of  the  liquid  in  the  groove  is  given  by 

Dt2 


Po  =  /Re 


2V 


(56) 


The  dimensionless  hydraulic  diameter  for  the  flow  of  liquid  in  a  sinusoidal  groove  with 
a  circular  meniscus  is  D ^  =  4 A^/P*,  where  the  dimensionless  cross-sectional  area  of 
the  liquid  is  given  by 
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The  dimensionless  wetted  perimeter  is  given  by  the  following  integral  equation. 
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The  mean  velocity  is  defined  as 
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See  Appendix  B  for  derivations  of  the  dimensionless  Poisson  equation,  dimensionless 
boundary  conditions,  and  dimensionless  liquid  cross-sectional  area  and  perimeter. 


2.3  Numerical  Model 

The  elliptic  Poisson  equation  given  by  eqn.  (48)  with  mixed  boundary  conditions 
[eqns.  (49),  (50)  and  (54)]  was  solved  using  Gauss-Seidel  iteration  with  a  central 
differencing  scheme  and  successive  over-relaxation  (Anderson  et  ah,  1984).  The  con¬ 
vergence  criteria  for  the  iterative  solution  was  set  to  e  =  10-10  for  each  case.  A  grid 
independence  check  was  made  in  which  the  number  of  grids  in  each  direction  was 
doubled.  When  the  value  for  the  Poiseuille  number  did  not  change  by  more  than  1%, 
grid  independence  was  considered  to  be  reached.  The  convergence  criteria  was  then 
reduced  by  an  order  of  magnitude  while  maintaining  the  highest  number  of  grids.  If 
the  Poiseuille  number  did  not  change  by  more  than  1%,  the  solution  was  considered 
to  be  independent  of  both  grid  size  and  e.  Otherwise,  a  grid  independence  check  was 
made  at  the  smaller  value  of  e  until  a  converged  solution  was  reached.  The  integral 
equation  for  the  wetted  perimeter  [eqn.  (59)]  was  integrated  numerically  since  no 
closed-form  solution  exists. 

The  numerical  model  was  tested  against  an  existing  solution  in  the  archival  liter¬ 
ature.  Shah  (1975)  determined  the  friction  factors  for  the  laminar  flow  within  ducts 
of  various  cross  sections  using  a  least-squares-matching  technique.  Table  3  shows 
the  comparison  of  the  Poiseuille  number  between  the  present  solution  and  that  given 
by  Shah  (1975)  for  laminar  flow  in  a  family  of  sinusoidal  ducts.  The  agreement  is 
excellent,  with  a  maximum  difference  of  1.1%. 
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Table  3:  Poiseuille  number  versus  sinusoidal  duct  aspect  ratio:  Comparison  of  the 
present  solution  with  that  given  by  Shah  (1975). 


fi 

Poiseuille  Nu 
Shah  (1975) 

mber,  Po 
Present 

1/4 

14.553 

14.479 

1/3 

14.022 

13.931 

1/2 

13.023 

12.935 

l/\/3 

12.630 

12.543 

2/3 

12.234 

12.148 

1 

11.207 

11.115 

2 

10.123 

10.061 

4 

9.743 

9.6373 

2.4  Results  and  Discussion 

A  numerical  study  has  been  completed  in  which  the  flow  of  liquid  in  a  sinusoidal 
groove  has  been  solved.  Figure  23  presents  contour  plots  of  the  dimensionless  liquid 
velocity.  The  maximum  liquid  velocity  increases  with  cocurrent  shear,  and  decreases 
with  countercurrent  shear,  as  expected.  For  countercurrent  vapor  flow,  a  portion  of 
the  liquid  flows  in  the  —y  direction,  which  is  opposite  to  the  direction  of  the  pressure 
gradient.  This  flow  reversal  shows  the  potential  of  the  vapor  shear  to  drive  the  mean 
velocity  of  the  liquid  to  zero,  or  to  completely  reverse  the  flow,  depending  on  the 
magnitude  of  the  pressure  gradient. 

Figure  24  and  Table  39  show  the  mean  velocity,  volumetric  flow  rate  and  Poiseuille 
number  versus  shear  stress  at  the  liquid-vapor  interface  for  several  values  of  the 
meniscus  contact  angle.  The  mean  velocity  increases  with  both  r*v  and  <f> ,  but  is  a 
relatively  weak  function  of  <f>.  For  a  given  value  of  meniscus  contact  angle,  the  mean 
velocity  is  linear  with  Tj*  due  to  an  overall  force  balance  on  the  liquid.  The  volumetric 
flow  rate  also  increases  with  rt*  and  4>,  but  is  a  much  stronger  function  of  <f>  due  to  the 
increase  in  the  cross-sectional  area  of  the  liquid.  The  Poiseuille  number  decreases  as 
t*v  increases.  For  cocurrent  vapor  flow  (t*v  >  0),  Po  decreases  steadily  with  t*v.  For 
countercurrent  flow,  Po  increases  dramatically  with  —  t*v  due  to  the  mean  velocity 
approaching  zero.  In  general,  Po  increases  with  0  due  to  the  increase  in  the  hydraulic 
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Figure  23:  Dimensionless  velocity  fields  for  laminar  flow  in  a  sinusoidal  groove 
0.5,  w*/2  =  0.25,  <j>  =  25°):  (a)  t*v  =  2.0  (cocurrent  flow);  ( b )  rfv  =  0.0;  (c)  = 

(countercurrent  flow). 


diameter  of  the  liquid. 

Figure  25  and  Table  40  present  the  mean  velocity,  volumetric  flow  rate  and 
Poiseuille  number  versus  the  groove  fill  ratio  for  several  values  of  the  groove  aspect 
ratio.  The  mean  velocity  increases  monotonically  with  area  ratio  and  groove  aspect 
ratio.  The  volumetric  flow  rate  also  increases  with  w*/w*  and  /?,  but  drops  off  rapidly 
for  w* /w*  <  0.5  due  to  the  decrease  in  flow  area.  The  Poiseuille  number  is  relatively 
constant  for  this  case.  For  /?  >  0.5,  Po  decreases  and  then  increases  with  area  ratio. 
For  j3  =  0.25,  Po  increases  over  the  range  of  w*/w*  examined. 

2.5  Semi- Analytical  Solution  for  v* 

As  seen  in  Fig.  24(a),  the  mean  velocity  is  a  linear  function  of  the  imposed  shear 
stress  at  the  liquid-vapor  interface.  Since  a  direct  numerical  simulation  of  the  liquid 
flow  field  for  a  number  of  values  of  the  shear  stress  is  computer  resource  intensive, 
it  is  appropriate  to  seek  a  semi-analytical  expression  for  v*.  Figure  26(a)  shows  the 
definition  of  the  parameters  involved,  where  the  mean  velocity  when  the  shear  stress 
is  zero  (uq)  is  given  by  the  numerical  solution.  The  value  for  the  liquid- vapor  shear  for 
which  the  mean  velocity  is  zero  {t*v  0)  is  given  by  the  following  force  balance  analysis. 
Figure  26(6)  shows  a  differential  element  of  the  liquid  in  the  groove.  A  force  balance 
between  the  pressure  drop  and  the  shear  forces  at  the  liquid-vapor  interface  and  at 
the  wall  results  in  the  following  relation. 

Py  Al  PyJfdyAl  -(-  TlvAlv  twAw  0  (®^) 

The  areas  over  which  the  shear  stresses  t1v  and  tT  act  are  Atv  =  Pivdy  and  Aw  -  Pdy. 
respectively.  Using  these  areas  and  nondimensionalizing  gives 

a;  +  r, •„*>,;  -  ^P*  =  0  (62) 

For  Poiseuille  flow  in  ducts  of  arbitrary  cross  section,  and  combined  Couette-Poiseuille 
flow  between  flat  plates,  the  shear  stress  at  the  wall  is  related  to  the  mean  velocity  of 
the  fluid  by  a  constant  (White,  1991).  Therefore,  in  the  present  analysis,  it  is  assumed 
that  this  also  holds  for  the  flow  of  liquid  in  a  sinusoidal  groove  with  an  imposed  shear 
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Figure  24:  Variation  of  the  flow  variables  with  shear  stress  at  the  liquid- vapor  interface 
for  various  values  of  meniscus  contact  angle  ((3  =  0.5,  w*/ 2  =  0.25,  P*  =  1.15245): 
(a)  Mean  velocity;  ( b )  Volumetric  flow  rate;  (c)  Poiseuille  number. 
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Figure  25:  Variation  of  the  flow  variables  with  groove  fill  ratio  for  various  values  of 
groove  aspect  ratio  (t*v  =  0.0,  <f>  =  0°):  (a)  Mean  velocity;  (6)  Volumetric  flow  rate; 
(c)  Poiseuille  number. 


Figure  26:  Semi-analyticalsolution  for  v*:  (o)  Definition  of  parameters;  (b)  Force 
balance  on  the  liquid  in  a  sinusoidal  groove. 
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stress  at  the  liquid- vapor  interface. 


(63) 


It  should  be  noted  that  the  constant  C\  is  probably  a  function  of  the  groove  geom¬ 
etry  and  meniscus  contact  angle.  However,  since  the  objective  of  this  analysis  is  to 
determine  the  liquid-vapor  shear  stress  when  the  mean  liquid  velocity  is  zero,  this 
functionality  is  unimportant.  The  perimeter  of  the  liquid-vapor  interface  is 

K  =  2K'  sin-1  (H)  (64) 


Substituting  these  relations  into  the  force  balance  equation  results  in  the  following 
expression  for  mean  velocity. 


v*  = 


1 

CXP* 


A*  +  2t;vR*  sin-1 


(65) 


The  mean  velocity  is  zero  when 


rlv, 0  = 


A* 


2  R*  sin  1 


(66) 


Figure  27(a)  shows  the  results  of  eqn.  (66).  The  numerical  results  shown  in 
Fig.  24(a)  were  extrapolated  to  determine  the  values  for  shear  stress  at  the  liquid- 
vapor  interface  when  v*  =  0.  Both  curves  indicate  that  t*v  0  increases  with  <j>  due 
to  the  increasing  depth  of  liquid  in  the  groove.  The  maximum  difference  between 
the  numerical  and  analytical  data  is  34%.  The  equation  for  the  mean  velocity  as  a 
function  of  the  shear  stress  is  given  by 


v'  =  V*/Vq  =  1  —  t' 


(67) 


where  r'  =  0.  The  semi-analytical  solution  for  the  normalized  mean  velocity 

is  shown  in  Fig.  27(b)  with  the  corresponding  numerical  data  presented  in  Fig.  24. 
Equation  (67)  predicts  94%  of  the  data  to  within  ±20%  over  the  range  of  meniscus 
contact  angle  examined. 
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Figure  27:  Comparison  of  the  semi-analytical  solution  with  numerical  data  (f3  = 
0.5,  w*/2  =  0.25):  (a)  Countercurrent  vapor  shear  stress  required  for  v*  =  0;  (6) 
Normalized  mean  velocity  versus  normalized  shear  stress  at  the  liquid-vapor  interface. 
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2.6  Conclusions 


Based  on  the  results  of  the  numerical  model  of  the  flow  of  liquid  in  a  sinusoidal 
groove,  the  following  conclusions  have  been  made:  For  a  given  meniscus  contact  angle, 
the  mean  velocity  was  linear  with  imposed  shear  stress  at  the  liquid-vapor  interface. 
The  volumetric  flow  rate  in  the  groove  was  negligible  for  groove  fill  ratios  of  less  than 
w* /w*  <  0.5.  The  Poiseuille  number  was  a  strong  function  of  the  countercurrent 
shear  stress.  A  semi-analytical  expression  was  provided  to  approximate  the  mean 
velocity  as  a  function  of  the  shear  stress  at  the  liquid-vapor  interface. 
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Nomenclature 


Ag 

At 

Aiv 

Ayj 

A* 

(X'f 

A 

b 

Bo 

d* 

D 

Dh 

D'k 

Fl 

F2 

f 

9 

h 

hfg 

hi 

h\ 

i 

k 

k 

La 

Lc 

Le 

Leff 

Lp 


cross-sectional  area  of  the  groove,  m2 
cross-sectional  area  of  the  liquid,  m2 
area  of  the  liquid- vapor  interface,  m2 
area  of  the  groove  wall,  m2 
AJh? 

radial  acceleration  vector,  m/s2 

acceleration  vector  at  any  point  in  the  helical  groove,  m/s2 

distance  from  the  liquid- vapor  interface  to  the  bottom  of  the  groove,  m 

Bond  number,  pgb2 /a 

parameter  defined  in  eqn.  (6) 

diffusion  coefficient 

hydraulic  diameter,  4 AJP,  m 

Dh/h 

mean  velocity  parameter,  hfv*,  m2 

volumetric  flow  rate  parameter,  hfV*,  m4 

friction  coefficient,  2f^/pv2 

acceleration  due  to  gravity,  m/s2 

groove  height,  m 

heat  of  vaporization,  J /kg 

height  of  the  liquid  in  the  groove  at  the  wall,  m 

ht/h 

unit  vector  in  the  x  direction 
unit  vector  in  the  z  direction 
curvature  vector 
adiabatic  length,  m 
condenser  length,  m 
evaporator  length,  m 

effective  heat  pipe  length,  Le/2  +  La  +  LJ 2,  m 
pitch  length,  m 
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u 

m 

n 

n* 

h 

Ng 

P 

P 

Plv 
p * 

K 

Po 

Qcap 

Qg 

Qt 

R 

Rc 
Rh 
Rv 
R * 

r; 

Re 

t 

TsaX 


V 

V 

Vl,  max 


total  heat  pipe  length,  m 
mass  flow  rate,  kg/s 

coordinate  normal  to  the  liquid-vapor  interface 
n/h 

unit  vector  normal  to  the  liquid- vapor  interface 
number  of  grooves 
pressure,  N/m2 
wetted  perimeter,  m 

perimeter  of  the  liquid-vapor  interface,  m 
P/h 

Piv/h 

Poiseuille  number,  / Re,  D*^/ 2v* 
capillary  limit  heat  transfer,  W 
heat  transfer  due  to  a  single  groove,  W 
total  heat  transported,  W 

radius  of  curvature  of  the  liquid-vapor  interface,  m 

capillary  radius,  m 

radius  of  the  helix,  m 

radius  of  the  heat  pipe  vapor  space,  m 

R/h 

dimensionless  radius  of  curvature  at  bifurcation 
Reynolds  number,  pvDh/fj, 
time,  s 

saturation  temperature,  K 
rr-direction  velocity,  m/s 
y-direction  velocity,  m/s 
mean  y-direction  velocity,  m/s 
maximum  mean  liquid  velocity,  m/s 
maximum  mean  vapor  velocity,  m/s 
pv/h2(-dp/dy) 
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V* 

Vo 

v' 

V 

y* 

w 

w* 

Wi 

wt 

*>  y,  Z 

X * 

Z 1 
z* 

p 

e 

e 


M 

P 

a 

7~lv 

Tlv 

rlvfi 

Tw 

T* 

1  w 

t' 

<fi 

00 

4>b 

LO 


dimensionless  mean  y-direction  velocity 
dimensionless  mean  y-direction  velocity  when  r*v  =  0 

V*/Vq 

volumetric  flow  rate,  vAi,  m3/s 
fiV/lh^-dp/dy)] 

width  of  the  bottom  of  the  groove,  m;  ^-direction  velocity,  m/ s 
w/h 

width  of  the  liquid  in  the  groove,  m 
wi/h 

coordinate  directions 
xjh 

distance  from  the  liquid-vapor  interface  to  the  groove  bottom,  m 
z/h 

groove  aspect  ratio,  w/2 h 
convergence  criteria 
groove  half-angle,  rad 
absolute  viscosity,  Pa-s 
density,  kg/m3 
surface  tension,  N/m 

shear  stress  at  the  liquid- vapor  interface,  N/m2 
Tiv/h(-dp/dy) 

shear  stress  at  the  liquid- vapor  interface  when  v*  =  0 
average  shear  stress  at  the  wall,  N/m2 
T^/h(—dp/dy) 

r*v/T*v,0 

meniscus  contact  angle,  rad 
minimum  meniscus  contact  angle,  rad 
meniscus  contact  angle  at  bifurcation,  rad 
over-relaxation  parameter 
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A  Derivation  of  Selected  Expressions  for  Trapezoidal  Grooves 


A.l  Determination  of  dimensionless  Poisson  equation 

The  governing  equations  for  fluid  flow  are  the  conservation  of  mass  and  conserva¬ 
tion  of  momentum  equations.  The  conservation  of  mass  equation  is  given  by  White 
(1991). 


dp 

dt 


+  divpV  =  0 


(68) 


The  conservation  of  momentum  equations  (in  scalar  form)  are  also  given  by  White 
(1991) 


Du  dp  d  (  du  ,-A  d 

p-rs  =  PS,  -  -5T  +  31  (  +  A  divV  )  + 


Dt 


dx  dx 


dy 


d 

+  dz 


(  du  dv' 
P\dy  +  dx 


dw 

"i  &  + 


du\ 

dz)m 


Dv  dp  d  [  dv 

PT)i  =  P9y  ~  dy  +  dx  [P  V&r  + 


du 

dy 


d  (  dv 

+  Tv  (2^  +  A  dlvV 


Dw 


dp  ,  d 


p^i  =  p9‘~d~z  +  T 


dw  du 

P[T  +  Tz 


d 

+  dz 


d_ 

dy 


dv  dw\ 

M Tz  +  T)\ 

dv  dw\ 

p{Tz+T  jj 


(69) 


where 


d  ( .  dw  - 

+s  War+AdlvV 


X  +  -p  —  0 

Expanding  eqn.  (68)  and  assuming  incompressible  flow  ( p  =  constant)  gives 

/  du  dv  dw' . 

p{T  +  T  +  ^]=0 


and  upon  simplification 


du  dv  dw  _ 
dx  +  dy  +  dz 


(70) 


(71) 


(72) 
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-4 

Solving  for  A  in  eqn.  (70),  substituting  it  into  eqn.  (69),  and  expanding  divV  results 


=  P9x  ~  7T-  + 


dp  d  (  du  2  \du  dv  dw 


2^  -  XM 


dx  dx  \  ^  dx  o  \_dx  dy  dz 


+  t. - b 


d  fdu  dv 
+  dy  ^  \dy  +  dx 


+  ¥z  " 


dw  du 
dx  dz 


dp  d  \  ( dv  du 


dy  dx  P  dy 


dv  2  f  du  dv  dw 


Ty  \2%  “  3P  |&  T  dy  T  dz \)  '  dz  r  \dz  '  dy 


+  ^-  +  -=-  + 


( dv  dw 

Hr  +  r 


-is  uj 

P-RI  =  P9z~  -TZ  + 


dw  +  du 


dz  dx  ^  \dx  dz  J I  dy  I  \dz  dy 


d  ‘ 

+  sr.  P 


dv  dw 
- f-  — “ 


v  .  dw  2  \du  dv  dw 

+&  [di  +  d-y  +  -T\) 

Assuming  constant  viscosity,  substituting  eqn.  (72)  into  eqn.  (73),  and  expanding 
the  total  derivatives,  the  conservation  of  momentum  equations  can  be  written  as 


du  du  du  du\  dp 

M+urx+vT+wTz)=pa’-Tx+>‘ 

( dv  dv  dv  dx A  dp 

ym+udi+vT+wrz)=P9’-T+,‘ 


d2u  d2u  d2u 

dx2  +  dy2  dz2 

d2v  d2v  d2v 


P9y  dy  +  p  \dx2  +  dy2  +  dz2 


dw  dw  dw  dw 
—  +  u—  +  v-z-  +  w— 
dt  dx  dy  dz 


dp 

ps*-Tz+p 


d2w  d2w  d2w 
dx2  dy2  dz2 


Equation  (74)  can  be  simplified  further  by  assuming  steady  flow  (d/ dt  =  0)  and  that 
no  body  forces  are  acting  on  the  fluid  (gx  =  gy  =  gz  =  0).  This  simplification  results 


(  du  du  du 
p[udi  +  vT  +  aTz 


( d2u  d2u  d2u 

+  P  V  dx2  +  dy2  +  dz2 


(  dv  dv  dv 

P[UT+VT+Wrz 

(  dw  dw  dw N 

p\uT  +  v^  +  wT, 


dp 

=  -jf+p 
dy 


dz+P 


( d2v  d2v  d2v 

\dx2  +  dy2  +  dz2 

d2w  d2w  d2w 

dx2  dy2  dz2 
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Finally,  by  assuming  fully  developed  flow  (u  =  w  =  0,  dv/dy  —  0)  eqns.  (72)  and 
(75)  reduce  to 

dv 

Continuity  :  =  0 

dy 


Momentum  : 


d2v  d2v  _  1  dp 
dx 2  +  dz 2  y  dy 


(76) 


dp  _  dp  _  o 
dx  dz 

Because  v  does  not  change  with  y,  dp/dy  is  constant,  and  the  momentum  equation 
reduces  to  the  Poisson  equation. 


d2  v  d2v  1  dp 
dx2  +  dz2  ydy 


(77) 


The  following  definitions  of  dimensionless  variables  can  be  used  to  write  the  Poisson 
equation  in  dimensionless  form 


v  = 


Solving  for  v,  x,  and  z  gives 


v 


yv  ,  _  x  *  _  z 

h2(— dp/dy)  X  h  h 


^■—v*h2  x  —  x*h  z  =  z*h 


ydy 


Substituting  these  relationships  into  eqn.  (77)  gives 


1  dp 
ydy 


d2  (  ^  [  --^Vh2 


+ 


1  dp 
ydy 


1  dp 
ydy 


(78) 


(79) 


(80) 


d{x*h)  '  d  (z*h) 

Equation  (80)  can  then  be  simplified  to  the  Poisson  equation  in  dimensionless  form. 

(81) 


dV  <9V 
+ 


-1 


dx*2  dz*2 

A. 2  Determination  of  dimensionless  boundary  conditions 


Figure  1  (a)  and  (b)  shows  liquid  flow  in  a  trapezoidal  groove  with  the  definitions 
of  the  dimensions  used  in  the  analysis.  The  no-slip  boundary  condition  is  applied  to 
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the  groove  walls  as  shown  below. 

0  <  x  <  w/2,  z  =  0 

u  =  0:  ^  w/2  <  x  <  w/2  +  htan9,  z  =  (x  —  w/2)  cot#  for  9>  0  (82) 


x  =  w/2,  0<z<h  for  6  =  0 

Using  the  definitions  of  the  following  dimensionless  parameters 

fj,v 


* 

v  = 


*  _  x  „  _  £  _  w_ 

h2(—dp/dy )  h  h  2  h 


(83) 


the  dimensional  no-slip  boundary  condition  can  be  written  in  dimensionless  form  as 
0<x*<P,  z*  =  0 


u*  =  0:  <  f3  <  x*  <  P  +  tan  9,  z*  =  (x*  —  0)  cot  9  for  9  >  0 

x*  =  P,  0  <  z*  <  1  for  9  =  0 

Along  the  centerline  of  the  groove,  the  symmetry  condition  is  used 

dv 


(84) 


dx 


=  0:  x  =  0,  0  <  z  <  (h  +  d)  —  R 


Using  the  dimensionless  parameters  from  eqn.  (83)  and  noting  that 

*■-? 

the  symmetry  condition  is  cast  in  the  form  of 
dv* 


dx* 


=  0  :  x*  =  0,  0  <  z*  <  (1  +  d*)  -  R* 


The  shear  stress  boundary  condition  is  applied  to  the  liquid  vapor  interface. 


r\ 

a— -  =Tiv  :  0  <  x  <  w/2  +  htan.9,  z  =  (h  +  d)  -  Vi?2  -  x2 
dn 

Using  the  dimensionless  variable  definitions  in  eqns.  (83)  and  (86)  along  with 

*  n 
n  =h 

the  shear  stress  boundary  condition  is  shown  below  in  dimensionless  form. 


=  r»*  :  0  <  x*  <  P  +  tan  9 ,  z*  =  (1  +  d*)  -  V R*2  —  x 
dn*  1  -- 


*2 


(85) 


(86) 


(87) 


(88) 


(89) 


(90) 
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A. 3  Equation  of  the  meniscus  for  9  >  0 

As  shown  in  Fig.  28(g),  the  equation  of  the  circle  which  is  symmetric  about  the 
^  axis  is  given  by 

x2  +  (z  -  z0)2  =  R 2  (91) 


where  z0  is  the  location  of  the  center  of  the  circle  in  the  z  direction. 


zq  =  h  +  d 


-  =  ^  +  htan0 
2  2 


Zo  =  h  +  y  R2  -  +  h  tan  0^ 

Therefore,  the  equation  of  the  liquid-vapor  interface  is 


(92) 


(93) 


A  relation  must  be  found  between  the  radius  of  curvature  and  the  meniscus  contact 
angle.  At  Point  2  shown  in  Fig.  28(6),  the  slope  of  the  liquid-vapor  interface  is  given 
by 


dz  _  (h  —  e) 
dx  h tan  9 

The  length  of  the  line  from  Point  2  to  Point  3  is 


(94) 


g  =  \/h2~+h2  tan2  9  =  hy/ 1  +  tan2  9 


(95) 


Therefore,  using  the  law  of  sines, 

hV  1  +  tan2  9  _  e 
sin  7  sin  4> 


where 


0  +  7  +  4>  =  tt 


(96) 


(97) 
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| -  w/2 - 1 

Figure  28:  Determination  of  the  equation  of  the  circular  liquid-vapor  interface. 
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Combining  the  above  two  equations  gives 


e  =  hy/  1  +  tan2  6 


siruj) 


Lsin(7r  —  9  —  4>) 
The  slope  of  the  liquid- vapor  interface  at  Point  2  is 


dz 

dx 


=  cot  6  —  y/l  +  cot2  6 


sin  (j) 


(98) 


(99) 


[sin(7r  —  0  -  (j>)  J 

Taking  the  derivative  of  the  liquid-vapor  interface  equation  [eqn.  (93)]  with  respect 
to  x  gives  the  slope  at  Point  2  in  terms  of  the  radius  of  curvature  of  the  liquid-vapor 
interface 

+  htmd'j 


dz 

dx 


(100) 


^J~R2  -  +  ht&n$Sj 

Equating  eqns.  (99)  and  (100)  gives  the  radius  of  curvature  in  terms  of  the  meniscus 


contact  angle. 


R  =  (!  +  /itanfl)  Jl+  jcot0  -  y/l  +  cot2  6 


A. 4  Equation  of  the  meniscus  for  0  =  0 


sin  (j> 


Lsin(7r  —  6  —  <j>) 


) 


-2 


for  0>O 

(101) 


For  the  case  when  the  groove  half-angle  is  zero  [Fig.  28(c)],  the  equation  of  the 
circular  arc  simplifies  to 


+  <  z 


h  + 


'*-(*)■ 


=  R2 


(102) 


dz 

dx 


The  slope  of  the  curve  at  Point  2  is  given  in  terms  of  the  meniscus  contact  angle. 

=  cot  4>  (103) 

Differentiating  eqn.  (102)  with  respect  to  x  gives  and  evaluating  at  Point  2  [x  =  w/2, 
z  =  h)  gives 

(!) 


dz 

dx 


»■-(?)’ 


(104) 
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Equating  eqns.  (103)  and  (104)  gives  the  meniscus  radius  in  terms  of  the  contact 
angle. 

R  =  (^)  \/l  +  tan 20  (105) 

A. 5  Point  of  meniscus  bifurcation  for  0  >  0 


The  radius  of  curvature  at  bifurcation  is  given  in  terms  of  the  geometry  of  the 
groove  as 

Rb=  2  I  h  +  —  ( —  +  h ^an ^ )  I  (10®) 


h+h^  +  ht& n9)  J 

The  equation  of  the  circular  arc  for  this  case  is  given  by 

x2  +  (z-  Rb )2  =  R\  (107) 

The  slope  at  Point  2  is  given  by  eqn.  (99).  Differentiating  eqn.  (107)  with  respect  to 


x  gives 


dz 

dx 


—  -H  h  t  an  0 
2 _ 

Rb  —  h 


(108) 


sin  & 


(109) 


Equating  eqns.  (99)  and  (108)  gives 

Rb  =  h  +  (^-  +  h  tan  #)  1  cot  6  -  \/l  +  cot2  6  ,  . 

\  2  /  (  [_sin(7r  -  6  -  <pb). 

The  meniscus  contact  angle  at  bifurcation  can  be  found  as  a  function  of  the  groove 

geometry  by  equating  eqns.  (106)  and  (109) 


(j)b  =  tan  1 


sin# 


cos  6  + 


2  sin  0  (p  +  tan  0) 

!  —  (/?  +  tan  0)2  J 


-l 


—  cos# 


(110) 


A. 6  Point  of  meniscus  bifurcation  for  #  =  0 


For  the  case  when  the  half-angle  is  zero,  the  lower  limit  of  recession  is  provided  as 

follows.  The  equation  of  the  circular  arc  is  given  by  eqn.  (107).  The  slope  at  Point 

2  is  shown  in  eqn.  (103).  The  slope  at  Point  2  can  also  be  found  by  differentiating 

eqn.  (107)  with  respect  to  x  and  evaluating  the  result  at  Point  2. 

w 


dz  _  2 

dx  2  R  —  h 


(111) 
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Equating  eqns.  (99)  and  (111)  gives 


w 


R  =  h  +  —  tan  <f> 


(112) 


A.7  Determination  of  the  hydraulic  diameter 


In  order  to  present  the  results  of  the  computations  in  terms  of  the  Reynolds 
number,  the  hydraulic  diameter  of  the  flow  field  must  be  found 

4  At 


Dh  = 


(113) 


where  A;  is  the  cross-sectional  area,  and  P  is  the  wetted  perimeter. 


P  =  w  +  2h\/l  +  tan2  6 


(114) 


The  area  of  the  liquid  is  found  by  subtracting  the  area  of  the  segment  of  the  circular 
arc  from  the  trapezoidal  groove  area.  The  area  of  the  trapezoid  is 


At  =  h(w  +  h  tan  6) 


The  area  of  the  circular  segment  is  given  by 

„-i  (  d 


where 


Ag  =  R2  cos"1  (  -  )  -  dVR2  -  d2 
H  A 


d 


(115) 


(116) 


(117) 


c  w  *  ,  n 
—  —  — I-  h tan# 
2  2 


(118) 


d  =  y  R2  -  +  h  tan  6^ 

d2  =  R2  -  +h  tan  0^ 


(119) 

(120) 


A.,  =  R2  cos-1  \l  1  -  -^2  (|  +  htan0)2  -  +  htan@)  R2  -  (|  +  htan0) 


(121) 
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The  cross-sectional  area  of  the  fluid  is  given  by 


1  / w 


Ai  =  h(w  +  hta.nO)  -  R2  cos  1  \  1  -  (  tt  +  htan#) 

V  V2  '  (122) 

+  +  h tan 9^j  yR2  -  +  h tan 0^ 

The  hydraulic  diameter  for  the  flow  of  liquid  in  a  trapezoidal  groove  with  a  circular 


meniscus  is 


Dh  =  4  {  h(w  +  h  tan  9)  —  R2  cos  1  \l  1-  ^  +  h  tan  0 


1  /u; 


R  v  2 


+J?  ^  +  htan^j  y  1  -  i^  +  htan# 


-\ 

>  +  2h\/l  +  tan2  0^ 


A. 8  Determination  of  the  unit  normal  vector  to  the  liquid- vapor  inter- 


The  unit  normal  vector  to  the  liquid-vapor  interface  is  related  to  the  curvature 


vector  by 


For  z  =  f(x),  the  curvature  vector  is 

k  = - Z—-,(-z',l) 

[1  +  (*')2] 

The  magnitude  of  the  curvature  vector  is  given  by 


\k\  = 


[l  +  (*f] 


^21 3/2 


The  functional  form  of  the  circular  arc  in  terms  of  z  =  f  (x)  is 


z  =  Zq  ±  \lKl  —  x2 


For  the  present  analysis,  the  lower  part  of  the  circular  arc  is  of  interest.  Therefore 


The  first  derivative  of  2  with  respect  to  x  is 


2'  =  xu 


(129) 


The  second  derivative  of  2  is 


z"  =  x*u-3'2  +  u~1/2 


(130) 


Substituting  the  above  relations  into  eqn.  (125)  gives  the  curvature  vector. 

k  =  u1'2  (1  +  A-1)'1  (-xtT1/2, 1)  (131) 


The  magnitude  of  the  curvature  vector  is 

|u“1/2  (i2«_1  +  1)  I 

11  “  (1  +  x2u ~'f2 

Since  u  =  R2  -  x2  >  0,  the  above  equation  reduces  to 

\k\=u-1'2  (*2u-1+l)'1/2 
The  unit  normal  vector  to  the  liquid-vapor  interface  is 


(132) 


(133) 


A. 9  Derivation  of  the  closed  form  solution  for  the  capillary  limit 


(134) 


The  equation  for  the  maximum  capillary  limit  which  includes  the  effects  of  inter¬ 
facial  shear  stress  and  groove  fill  ratio  is  shown  below  for  one  helical  groove. 


o  >Teff 

Rc  hfg 


A/ivQt 


&PvQt  ,  f  piQ^ _ _ 

'  7 rpvRi  \hfv^piAi  7r hiRlpv 


-Pi 


(135) 


Assuming  no  body  forces,  the  previous  equation  reduces  to 

4  (J'vQt 


&  >  ^eff 


Rr 


lfg 


( 


2r;Rh 


SpvQt  /  PdQg _ 

7 rpvR*  yhfvQpiAi  irhiR^Pv  rtv,oJ  V  V  LP 


+  1 


(136) 
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By  allowing  Lp  -*  oo,  which  assumes  the  grooves  are  straight  axial  grooves  rather 
than  helical,  eqn.  (136)  can  be  written  as 

°  >^eff  ( WQg _ fyjLyQt 

Rc  -  hfg  yrpvRi  hfv^piAt  TrhiRlpvT^v0 

Factoring  terms  on  the  right  hand  side  of  the  equation  gives 

<7  ^  Lefi  P’lQg  &P>vQt  |  2  _  R^_ 

Rc  hfg  hfvoPiAi  TrpvRy  y  2htTlvfi 

The  total  heat  transfered  by  the  heat  pipe  is 

Qt  =  'E<i,  (139) 

2=1 

Assuming  each  groove  transports  the  same  amount  of  heat,  eqn.  (139)  can  be  sim¬ 
plified  to 

Q,  =  f)  (140) 

Substituting  this  relationship  into  eqn.  (138)  gives 

<7  ^I'eff  PdQt  |  8 PvQt  Ry 

Rc  “  hfg  NghfolpiAi  KpvRt  y  2h;rz*  0 

Solving  for  <3*,  the  resulting  equation  is  the  maximum  capillary  limit  in  a  heat  pipe 
acting  with  no  body  forces  and  straight  axial  grooves  with  the  effects  of  interfacial 
shear  stress  and  groove  fill  ratio. 

•  _  ahfg  Pi  &Pv  /,  Ry  ^  (142) 

^cap  "  RcLrt  [NghffiptAi  +  tt pvRi  ^  2 htr^  ) 
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B  Derivations  of  Selected  Expressions  for  Sinusoidal  Grooves 


B.l  Determination  of  dimensionless  Poisson  equation 


The  governing  equations  of  fully  developed  laminar  flow  reduce  to  the  Poisson 
equation,  as  shown  by  White  1991. 


d2v  d2v  1  dp 

+ 


(143) 


dx 2  dz 2  p  dy 

By  substituting  the  following  definitions  for  the  dimensionless  variables  used  in  the 
analysis 

pv 


______  *  £  *  _  £ 

V  h2{—dp/dy)  h  h 


(144) 


eqn.  (143)  can  be  rewritten  as 


d2v*  d2v* 


=  -1 


dx*2  dz*2 

B.2  Determination  of  dimensionless  boundary  conditions 


(145) 


Figure  22  shows  liquid  flow  in  a  sinusoidal  groove  with  the  definitions  of  the 
dimensions  used  in  the  analysis.  The  no  slip  boundary  condition  is  applied  to  the 
groove  walls  as  shown  below. 

2x 


=Mi+cosbS+l)]} 


v  =  0  :  0  <  x  <  wij 2,  z 

Using  the  definitions  of  the  following  dimensionless  parameters, 

„  pv  *_£  *  _  £  o=i£ 

V  h2(—dp/dy )  h  h  2  h 


(146) 


(147) 


the  dimensional  no  slip  boundary  condition  can  be  written  in  dimensionless  form  as 


v*  =  0  :  0  <  x*  <  w*/2,  z*  =  -  <|  1  +  cos 


x’ 


"7rf7  +  1 


)1) 


(148) 


Along  the  centerline  of  the  groove,  the  symmetry  condition  is  used 


dv 

dx 


=  0  :  x  =  0,  0  <  z  <  hi  +  R^J  1  —  (^)  “  R 


(149) 
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Using  the  dimensionless  parameters  from  eqn.  (147)  and  noting  that 

*•-! 

the  symmetry  condition  is  cast  in  the  form  of 


(150) 


'  |)  =  0:  x'  =  0,  0<z'<hl  +  R^l-(j^j  _*•  (] 

The  shear  stress  boundary  condition  is  applied  to  the  liquid  vapor  interface. 

=  rlv  :  0<x<  wif  2,  z  —  hi  +  RyJ 1-  -  VR2-x2  (1 

Using  the  dimensionless  variable  definitions  in  eqns.  (147)  and  (150)  along  with 


"  =z 


(153) 


the  shear  stress  boundary  condition  is  shown  below  in  dimensionless  form. 


^  =  r;v:  0<x*<wU2,z*=h:  +  R^l-(j^j  -  V-R*2  -  z*2  (154) 

B.3  Equation  of  the  meniscus 


To  determine  the  equation  of  the  circular  meniscus,  a  relationship  must  be  made 
based  on  the  slope  of  the  tangent  line  at  the  sinusoidal  groove  wall,  which  is  given  by 

tan  ijj  =  (155) 

dx*  sine 

The  equation  of  the  sinusoid  for  the  groove  wall  is  given  by 

2*  =  £  jl  +  COS  -7T  +  1^  |  (156) 

Differentiating  eqn.  (156)  with  respect  to  x*  gives 

5?  =  ^sin[-,r(j  +  1)]  (157) 

Therefore,  substituting  eqn.  (157)  into  eqn.  (155)  and  evaluating  x*  =  w*/ 2  results 
in  the  equation  of  the  slope  of  the  tangent  line  at  the  wall  of  the  sinusoidal  groove. 


tam/> 


7T  f  WT 
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R* 


Figure  29:  Definition  of  geometric  parameters. 


The  slope 


From  Fig. 


of  the  tangent 


29,  the  angle 


line  on  the  circular  meniscus  is  given  by 
dz* 


tan  9  = 


dx* 


circle 


9  can  be  found  as  a  function  of  angles  a,  <f>,  and  ip 


a  =  7 r  —  ip 
6  +  (p  +  a  =  7r 
9  =  ip  —  <p 


(159) 


(160) 

(161) 

(162) 


Taking  the  tangent  of  both  sides  of  eqn.  (162)  and  using  the  trigonometric  addition 
formula  for  tangent  gives 

tan  9  =  tan  (ip  -  (p)  (163) 


tan(a  ±  /?) 


tan  a  ±  tan  (3 
1  tan  a  tan  (3 


(164) 
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(165) 


.  tan  ip  -  tan  6 

tan  9  =  - - - - T 

1  +  tan  ip  tan  <p 

The  general  equation  for  a  circle  with  the  z*  axis  as  a  line  of  symmetry  is 


x*2  +  (z*  -z*0)2  =  R *2 

Taking  the  derivative  of  eqn.  (166)  with  respect  to  x*  gives 

dz* 

2l*  +  2 (*•  -  z-)—  =  0 

Rearranging  eqn.  (167)  to  solve  for  dz’  /dx*  gives 


(166) 


(167) 


T~.  =  -7-T— 7 7  <168> 

dx*  {z*  —  z*) 

By  substituting  x*  =  w*/2  and  z*  =  h*  into  eqn.  (168)  and  then  substituting  the 

result  into  eqn.  (159)  gives  the  slope  of  the  tangent  line  on  the  circular  meniscus 

dz*  w*/2  ,  Qs 

TT  =  = tan0  (169) 

dx  circle  ~  Zo) 

Combining  eqns.  (165)  and  (169)  gives 


(169) 


Wi/2  tan^/i  —  tan<^ 
(h*  —  z*)  1  +  tan  ip  tan  <p 


(170) 


Rearranging  eqn.  (170)  to  solve  for  z*  gives 


z*0  =  h*l  + 


(1  +  tan  ip  tan  <p) 
tan  ip  —  tan  </> 


where  z*0  is  the  height  of  the  center  of  the  circle  on  the  2*  axis.  By  letting  d  be  the 
distance  between  the  center  of  the  circular  liquid-vapor  interface  and  the  height  of 
the  liquid,  it  can  be  shown  that 


z*0=h\  +  d 


(172) 


(173) 


2  2 
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(175) 


<  =  '»:+ 


2 


Combining  eqns.  (171)  and  (175)  gives 


J  (1  +  tan'i/'tan^) 
tan  xp  —  tan  <j> 


R*2 


^  (1  +  tan  ?/>  tan  </>) 
tan  xp  —  tan  <p 


(176) 


(177) 


where 


(178) 


7T 


d*  =  tan  ip  =  —  sin 
2/5 


-7T 


^  +  1 

2)0 


(179) 


/  Oil*  \ 

I  -r- )  (1  +  d*  taxup) 

\  £  J 

d *  —  tan  (p 


(180) 


R*2  - 


/  *\  2 

(  J  (1  +  d*  tan  (p)2 
(d*  —  tan0)2 


(181) 


i? 


»2 


(1  +  d*  tan  (p)2 
+  (d*  —  tan  (py ) 


(182) 


The  equation  of  the  radius  of  curvature  of  the  meniscus  in  terms  of  meniscus  contact 


angle  is 


(1  +  d*tan0)2"  2 
_  +  (d*  -  tan  (p )2) 


(183) 
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Figure  30:  Determination  of  the  cross-sectional  area  of  the  liquid. 

B.4  Determination  of  the  maximum  meniscus  contact  angle 

The  maximum  meniscus  contact  angle  can  be  found  when  R*  oo  in  eqn.  (183). 
It  can  be  seen  that  R*  -4  oo  as  (d*  —  tan  <f>)  — >  0.  Therefore 

d*  -  tan  0max  =  0  (184) 

Substituting  eqn.  (184)  into  eqn.  (179)  for  d*  results  in 

sin  -7T  +  1^  -  tan  <£max  =  0  (185) 

Solving  for  <frmax  gives  the  equation  for  the  maximum  meniscus  contact  angle. 

=  tan-'  {  Yp sin  [-*  (i  + 1)]  }  (186) 

B.5  Determination  of  the  hydraulic  diameter 

The  hydraulic  diameter  of  the  flow  field  must  be  determined  so  that  the  Poiseuille 
number  can  be  calculated.  To  find  D*h  =  4 A*/P*,  it  can  be  shown  from  Fig.  30  that 

Ai  +  A2  +  A$  =  —  w*h*  (187) 
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(188) 


*>  =  §  =  -  M  -  4. 

The  area  A3  can  be  found  by  integrating  the  equation  for  the  sinusoidal  groove  wall, 
eqn.  (156)  over  0  <  x*  <  w*/2. 


^3 


-L 


■W’J2  y  ( 

-  <  1  +  cos 


■*  (j+1)]}dI‘ 


(189) 


A  W*  P  . 

A3  =  T~2iSm 


wT 

~*'Tp  +  l 


(190) 


From  Fig.  30,  A\  is  the  area  above  the  meniscus.  It  can  be  shown  that  A\  is  one-half 
of  the  area  of  the  segment 

d 


A  R*2 
Ax  =  —  cos 


-l 


R' 


dVR*2  -  d2 


(191) 


where 


<!=,/**’- (f- 


(192) 


Thus,  the  area  A\  is 


,  _R* 


R*  cos 


-l 


'1- 


wj_ 
2  R* 


Wj 

2 


~~  I  A  /  1 


wl 

2R’ 


(193) 


Substituting  eqns.  (190)  and  (193)  into  eqn.  (188)  gives  the  cross-sectional  area  of 
the  flow  field 


wl 


A\=-^{2h\-l)-R* 


T)t-\2R 


w 


^  sin  -it  ^7^  +  1^  for  p  < 


(194) 


When  the  meniscus  contact  angle  approaches  the  maximum  value  <j)  ->  pma.x,  the  area 
Ai  — >■  0.  Therefore,  the  cross-sectional  area  of  the  flow  field  for  this  case  is 


A*  =  w^h]  -  2A3 


(195) 


A*  =  w*  hi  -  -  sin 

‘  ‘  2  7 T 


.  wl 

-%{W  +  1 


(196) 
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for  (j)  —  (j)  max 


(197) 


4’  =  ^  (2ftf-l)  +  §sin 

Z  7T 


W? 

-w'w+1 


In  order  to  determine  the  wetted  perimeter,  the  following  formula  will  be  used  to  find 
the  arc  length  of  the  groove  wall 


rx  2  _ 

P*  =  y/l  +  Jjffidx  for  y  =  f(x ) 
J  X 1 


where 


/•wf/2  / 

=2l  V1+ 


dz* 

dx* 


dx* 


7T  . 

(x*  ,  \1 

=  —  sm 

—  7T 

sine 

KP  )\ 

(198) 


dz :* 
dx* 

Thus,  the  arc  length  can  be  found  by  substituting  eqn.  (199)  into  eqn.  (198) 

2 


(199) 


=  2f 


7T 


'1  +  1  55'  sin 


-lj+1 


da:* 


(200) 


The  equation  for  the  wetted  perimeter  has  no  closed-form  solution,  and  consequently 
the  equation  for  the  hydraulic  diameter  has  no  closed-form  solution.  The  reason  is 
that  the  integral  in  eqn.  (200)  is  elliptic,  and  must  be  solved  numerically. 
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J  ~lj  1 
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1 - Ax* 

*  vV 1 

Figure  31:  Computational  molecule  using  a  central  differencing  scheme. 


C  Numerical  Code  Validation  (Trapezoidal  Grooves) 

The  elliptic  Poisson  equation  given  in  eqn.  (3)  with  mixed  boundary  conditions 
[eqns.  (4),  (5)  and  (8)]  was  solved  using  Gauss-Seidel  iteration  with  successive  over¬ 
relaxation  and  a  central  differencing  scheme  (Burden  and  Faires,  1985).  The  compu¬ 
tational  molecule  is  shown  in  Fig.  31.  The  second-order  central  differences  in  the 
x*  and  z*  directions  are 

av 


_ _  vi+u  Z  +  vJzlA 

dx*2  ~  (Ax*)2 

dV  ^  Vjj+i  ~  +  Vjj-i 

dz*2  ~  (A z*)2 


(201) 


(202) 


Substituting  eqns.  (201)  and  (202)  into  eqn.  (3)  gives  the  following  finite  difference 
equation 


*(fc) 

y.\  J  = 


+ 


|_(Ax*)2  (Az*)2J 


-1 


1  />(*-!>  , 
x*)2  +  J 


r  '  (Ax*) 

.  I- 

z*)2  vl- 


+  v*{k) 

+  vi,j- 1 


(203) 


(Az*)2  ViJ+l 

where  (k  -  1)  signifies  the  solution  field  from  the  previous  iteration.  Successive  over¬ 
relaxation  was  implemented  using  the  following  relation 

(204) 
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where  1  <  u)  <  2.  For  the  present  problem,  a  value  of  co  =  1.3  optimized  the 
convergence  of  the  solution.  Due  to  the  nature  of  the  rectangular  grid  used  to  solve 
eqn.  (203),  the  shear  stress  boundary  condition  at  the  liquid- vapor  interface  given  by 
eqn.  (8)  was  approximated  as  follows. 


dv*  dv*  * 
fa*  ~  fa*  ~Tlv 


(205) 


Discretizing  this  equation  gives  the  dimensionless  velocity  at  the  liquid-vapor  inter¬ 
face. 


*(fc)  *(&)  |  *  a 

=  vi,j- 1  +  Tiv 


(206) 


The  solution  procedure  was  as  follows: 

1.  Set  all  parameters. 

2.  Compute  the  radius  of  curvature,  area  and  hydraulic  diameter  of  the  flow  field. 

3.  Calculate  the  grid  sizes  and  locate  the  grid  points. 

4.  Locate  the  circular  liquid-vapor  interface  and  the  side  wall  of  the  groove. 

5.  Set  the  initial  velocity  field  to  zero. 

6.  Begin  iterating. 

7.  Calculate  the  updated  velocity  field  away  from  the  boundaries  using  eqns.  (203) 
and  (204). 

8.  Apply  the  symmetry  boundary  condition  at  x*  =  0. 

9.  Apply  the  shear  stress  condition  at  the  liquid- vapor  interface  using  eqn.  (206). 

10.  Set  the  velocity  to  zero  above  the  meniscus. 

11.  Calculate  the  difference  between  the  previous  solution  and  the  most  recent 
solution  at  every  grid  point. 
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(a)  If  the  maximum  difference  is  greater  than  the  convergence  criterion  e,  then 
place  the  most  recent  values  into  the  previous  solution  and  continue  to 
iterate. 

(b)  If  the  maximum  difference  is  less  than  the  convergence  criterion,  then 
calculate  the  average  velocity,  Poiseuille  number,  and  volumetric  flow  rate 
and  print  the  results. 

The  convergence  criteria  for  the  iterative  solution  was  set  to  e  =  10~8  for  each  case.  A 
grid  independence  check  was  made  in  which  the  number  of  grids  in  each  direction  was 
doubled.  When  the  value  for  the  Poiseuille  number  did  not  change  by  more  than  3%, 
grid  independence  was  considered  to  be  reached.  The  convergence  criteria  was  then 
reduced  by  an  order  of  magnitude  while  maintaining  the  highest  number  of  grids.  If 
the  Poiseuille  number  did  not  change  by  more  than  2%,  the  solution  was  considered 
to  be  independent  of  both  grid  size  and  e.  Otherwise,  a  grid  independence  check  was 
made  at  the  smaller  value  of  e  until  a  converged  solution  was  reached.  In  all  of  the 
results  reported,  grid  independence  and  the  convergence  criteria  were  less  than  1% 
unless  otherwise  noted. 

The  numerical  model  was  tested  against  several  existing  solutions,  such  as  rectan¬ 
gular  ducts,  triangular  grooves  without  interfacial  shear  stress,  and  rectangular  and 
triangular  grooves  with  shear  stress  at  the  liquid-vapor  interface. 

Shah  (1974)  determined  the  friction  factors  for  the  laminar  flow  within  ducts 
of  various  cross  sections  using  a  least-squares-matching  technique.  For  trapezoidal 
ducts,  the  area  and  hydraulic  diameter  are  given  by 

A*  =  2(3  +  tan  9  (207) 

=  2(2/?  +  tanf?)  (208) 

h  2  0  +  tan  9  +  sec  9 

Figure  32  and  Table  4  show  the  comparison  of  the  Poiseuille  number  between  the 
present  solution  and  those  given  by  Shah  (1975)  and  Shah  and  London  (1978)  for 
laminar  flow  in  a  family  of  trapezoidal  ducts.  The  agreement  is  excellent,  with  a 
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Figure  32:  Po  versus  (3  for  laminar  flow  in  a  rectangular  duct. 


maximum  difference  of  0.9%.  While  Shah  provided  a  Poiseuille  number  at  f3  =  0,  this 
could  not  be  obtained  with  the  present  numerical  solution  technique  because  the  area 
of  the  computational  domain  approached  zero  as  (3  — >  0.  As  a  result,  a  data  point  is 


shown  at  f3  =  0.01. 

Ayyaswamy  et  al.  (1974)  presented  the  friction  factors  obtained  for  laminar  flow 
in  triangular  grooves  using  the  Galerkin  method  of  solution.  Romero  and  Yost  (1996) 
presented  the  following  equation  for  the  dimensionless  volumetric  flow  rate  in  a  tri¬ 
angular  groove 


r(*,aO«r(a,aO  (h*c(<j>,a)Y 


3  ( A*((t>,a)\1/2 


cot  a 


(209) 


where  a  =  vr/2  -  6  is  the  angle  from  the  side  of  the  groove  to  horizontal,  and  r(a,  a) 
was  approximated  using  asymptotic  methods  and  a  regression  analysis. 

cot3  a  +  3.4  cot4  a  +  cot5  a 


T{a 


(210) 


(211) 


+  3.4  cot  a  -I-  4  cot2  a  +  3.4  cot3  a  +  cot4  a 

The  dimensionless  height  of  the  liquid  at  the  line  of  symmetry  of  the  groove  is 

_  cot  a  cos(o!  -  <f>)  -  1 
c  sin(o!  —  <j>) 

The  dimensionless  cross-sectional  area  of  the  liquid  in  the  groove  is 

^  2  [sin2 (ct  -  0)  tan  a  -  (ex  -  0)  +  sin(a  -  </>)  cos(o;  -  </>)] 

A  —  tan2 a  sin2 (a  -  4>) 
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Table  4:  Poiseuille  number  versus  rectangular  duct  aspect  ratio:  Comparison  of  the 
present  solution  with  that  given  by  Shah  (1975)  and  Shah  and  London  (1978). 


Poiseuille  Number,  Po 

0 

Shah  (1975) 

Shah  and  London  (1978) 

Present 

0.0 

24.000 

0.01 

23.619 

0.125 

20.585 

20.486 

0.2 

19.07050 

18.958 

0.25 

18.233 

18.116 

0.4 

16.36810 

16.246 

15.548 

15.424 

MUM 

14.97996 

14.856 

0.75 

14.476 

14.352 

0.8 

14.37780 

14.255 

1.0 

14.227 

14.22708 

14.104 

For  comparison  with  the  present  results,  the  Poiseuille  number  in  terms  of  the  dimen¬ 
sionless  volumetric  flow  rate  is 


d*2A* 
Po  =  h 


2m  a)  (213) 

Kolodziej  et  al.  (1999)  used  the  boundary  collocation  method  to  solve  the  same 
problem,  except  that  the  liquid- vapor  interface  was  not  assumed  to  be  circular.  For 
comparison,  the  case  where  Bo  =  0.001  was  examined.  Figure  33  and  Table  5  present 
the  solutions  obtained  by  Ayyaswamy  et  al.  for  9  =  5°  and  60°  for  the  full  range  of 
meniscus  contact  angle  (0.1°  <  <t>  <  9  +  4>  =  90°).  While  the  agreement  is  excellent 
for  9  =  60°,  it  is  less  so  for  9  =  5°.  This  is  due  to  the  extreme  narrowness  of  the 
flow  field  for  this  case.  The  maximum  percent  differences  for  9  =  5°  and  60°  are  3.7% 
and  0.9%,  respectively.  Figure  33  and  Table  5  also  show  a  comparison  of  the  present 
results  with  those  of  Romero  and  Yost  (1996).  The  maximum  percent  differences  are 
2.2%  for  9  =  5°  and  2.3%  for  9  =  60°.  Also  shown  in  Fig.  33  and  Table  5  are  the 
results  of  Kolodziej  et  al.  for  the  lowest  value  of  Bond  number  presented  (Bo  =  0.001), 
which  should  correspond  to  the  present  case.  The  agreement  with  the  Kolodziej  et 
al.  data  is  less  satisfactory,  with  a  maximum  percent  difference  of  4.5%  for  9  =  5°, 
and  19.9%  for  9  =  60°.  This  may  be  due  to  the  approximate  nature  of  the  solution 
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Figure  33:  Po  versus  <j>  for  laminar  flow  in  a  triangular  groove  without  shear  at  the 
meniscus  (/3  --  0,  t*v  =  0). 

by  Kolodziej  et  al.,  which  was  in  terms  of  a  truncated  infinite  summation. 

DiCola  (1968)  presented  the  solution  for  the  Poiseuille  number  for  the  laminar 
flow  of  a  constant  property  fluid  within  a  rectangular  groove.  While  interfacial  shear 
stress  at  the  liquid-vapor  interface  was  accounted  for,  the  groove  was  assumed  to  be 
completely  full,  with  a  meniscus  contact  angle  of  <f>  =  90°.  The  analytical  solution 
was  obtained  using  separation  of  variables.  Schneider  and  DeVos  (1980)  presented 
the  following  equation  which  was  attributed  to  DiCola  (1968): 

(2n  +  1)  7 t5 


24  <1  + 


(/Re),  = 


T*W2 
12 ¥ 


>-Sx: 


+ 1) 


sech- 


w 


/  w \2 

(1+2{) 


1 


96w 


E 


tanh 


(2 n  +  1)  tt5 


w 


(214) 


n$S  (2n  +  l)5 

where  the  dimensionless  shear  stress  at  the  liquid-vapor  interface  is  given  by 

Srv 


t*  = 


VlVl 


(215) 


The  groove  width  and  depth  are  w  and  6,  respectively,  tv  is  the  shear  stress  on  the 
vapor  at  the  liquid-vapor  interface,  [M  is  the  absolute  viscosity  of  the  liquid,  and  vi  is 
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Table  5:  Poiseuille  number  versus  meniscus  contact  angle  for  triangular  grooves  with 
no  shear  at  the  liquid- vapor  interface  (0  =  0,  t*v  =  0):  Comparison  of  the  present 
solution  with  those  given  by  Ayyaswamy  et  al.  (1974),  Romero  and  Yost  (1996)  and 
Kolodziej  et  al.  (1999). 


Poiseuille  Number,  Po 

Ayyaswamy  et  al. 

Romero  and  Yost 

Kolodziej  et  al. 

Present 

<t> 

II 

C n 
o 

_ 

0  =  60° 

0  =  5° 

II 

Ob 

o 

o 

0  =  5° 

0  =  60° 

o 

lO 

II 

o 

O 

CO 

II 

0.1° 

11.422 

7.612 

11.234 

7.645 

11.008 

7.545 

0.5 

11.436 

7.767 

11.248 

7.826 

1 

11.452 

7.963 

11.265 

8.046 

IKilil 

5 

9.403 

wnmm 

9.586 

10 

feWYil 

10.803 

HUH 

11.040 

10.803 

9.158 

11.291 

10.789 

15 

11.801 

11.841 

11.618 

12.093 

10.904 

10.215 

11.366 

11.812 

20 

11.884 

12.642 

11.704 

11.446 

12.606 

25 

11.952 

13.288 

11.777 

11.059 

N/A 

11.519 

13.397 

30 

12.007 

13.828 

emu 

11.120 

12.352 

11.581 

13.770 

35 

12.052 

11.171 

11.641 

KB 

11.931 

11.214 

45 

12.117 

■wrfiif 

50 

12.139 

11.999 

IB 

12.026 

IUB 

12.050 

IB 

65 

12.171 

12.071 

11.929 

— 

utilrll 

12.090 

11.969 

mm 

Emm 

12.108 

12.009 

80 

12.156 

12.124 

12.053 

85 

12.140 

12.063 
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Figure  34:  Po  versus  /?  for  laminar  flow  in  a  full  rectangular  groove  with  shear  at  the 
meniscus  (<f>  =  90°,  9  =  0°). 


the  mean  liquid  velocity.  The  DiCola  equation  can  be  rewritten  as 


Po  = 


J,_  tJ£  fl  96  V  1  :cch(2"+1),r 

1  3tr  L  ir>  So  (2n+i)4  2(3 

(1  +  BY  1 - i-  >  - r  tanh  1 — — TL— 

t  £s(2n+1)  2/5 


(216) 


A  negative  sign  has  been  introduced  to  account  for  the  fact  that  the  shear  stress  on 
the  vapor  is  in  the  opposite  direction  compared  to  that  on  the  liquid.  The  comparison 
between  the  above  equation  and  the  results  of  the  present  analysis  are  shown  in  Fig. 
34  and  Table  6  for  laminar  flow  in  a  family  of  rectangular  grooves  at  the  full  groove 
condition  (<£  =  90°,  0.1  <  0  <  1.0,  t{v  =  -0.1,  0.0,  and  1.0).  The  comparison  is 
excellent  with  a  maximum  percent  difference  of  2.3%. 

The  results  of  the  present  model  were  compared  to  the  correlation  presented  by 
Lin  and  Faghri  (1997),  where  the  friction  factor  for  the  flow  of  liquid  in  triangular 
grooves  was  presented  as  follows: 


fi  = 


13.163 

Re; 


+  0.6211t* 


(217) 


where  r*  =  Tv/pivi2  is  the  dimensionless  interfacial  shear  stress.  The  prediction  errors 
for  this  equation  were  between  -2.8%  and  1.3%.  Casting  this  equation  in  terms  of 
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Table  6:  Poiseuille  number  versus  rectangular  groove  aspect  ratio  for  several  values  of 
dimensionless  shear  stress  at  the  liquid- vapor  interface  (cf>  —  90°,  6  =  0°):  Comparison 
of  the  present  solution  with  that  given  by  DiCola  (1968). 


Poiseuille  Number,  Po 

p 

1—1 

O 

o 

II 

*  ^ 

)iCola  (196 

Tlv  =  1-0 

>8) 

Tlv  =  -0-1 

r,;  =  0.0 

Present 
^  =  10 

V,  =  -o.i 

0.1 

21.168 

10.300 

23.708 

21.199 

10.183 

23.788 

0.2 

19.070 

8.9579 

21.535 

18.952 

8.8446 

21.557 

0.3 

17.512 

7.9693 

19.945 

17.309 

7.8205 

19.940 

0.4 

16.368 

7.2283 

18.776 

16.219 

7.1135 

18.740 

0.5 

15.548 

6.7168 

17.934 

15.460 

6.6213 

17.891 

0.6 

14.979 

6.3590 

17.349 

14.874 

6.2958 

17.303 

0.7 

14.605 

6.1419 

16.966 

14.485 

6.0661 

16.920 

0.8 

14.377 

5.9964 

16.738 

14.245 

5.9265 

16.693 

0.9 

14.260 

5.9513 

16.630 

14.117 

5.8172 

16.586 

1.0 

14.227 

5.9070 

16.612 

14.158 

5.7804 

16.570 

the  dimensionless  quantities  used  in  the  present  analysis  gives 

Po  =  13.163  +  0.6211  [^r1 

Figure  35  and  Table  7  show  the  comparison  of  the  regression  equation  and  the  results 
of  the  present  analysis  for  countercurrent  flow.  Unfortunately,  not  enough 

information  was  provided  by  Lin  and  Faghri  (1997)  to  precisely  determine  the  limits 
of  applicability  for  eqn.  (217).  Therefore,  the  equation  was  evaluated  over  a  fairly 
large  range  for  comparison  with  the  present  solution.  The  agreement  is  quite  good 
between  0.075  <  -r*v  <  0.1,  where  the  maximum  percent  difference  in  this  range  is 
2.2%. 
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Figure  35:  Po  versus  —t*v  for  laminar  countercurrent  flow  in  a  triangular  groove 
(ft  =  0,  4>  =  30°,  9  =  20°). 


Table  7:  Poiseuille  number  versus  dimensionless  shear  stress  at  the  liquid- vapor  in¬ 
terface  for  countercurrent  flow  ((3  =  0,  (/)  =  30°,  6  =  20°) :  Comparison  of  the  present 
solution  with  that  given  by  Lin  and  Faghri  (1997). 


Poiseuille  Number,  Po 

msm 

Lin  and  Faghri  (1997) 

Present 

wmm 

13.862 

12.607 

0.050 

mMMM 

0.075 

15.733 

irarffip 

0.100 

17.026 

17.415 

0.125 

18.716 

0.150 

m 

0.175 

23.850 

27.533 
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D  Tables  of  Results 


In  all  of  the  results  reported,  grid  independence  and  the  convergence  criteria  were 
less  than  1%  unless  otherwise  noted. 


Table  8:  Mean  velocity  versus  shear  stress  at  the  liquid-vapor  interface  for  various 
values  of  meniscus  contact  angle  (j3  =  1.0,  0  =  0°). 


V* 

II 

O 

o 

o 

O 

CO 

II 

<t>  =  60° 

o 

O 

Ob 

II 

-e- 

— 

— 

— 

0.49199E-1 

— 

— 

— 

0.59164E-1 

— 

— 

0.32345E-1 

0.69436E-1 

— 

— 

0.42262E-1 

0.79708E-1 

— 

— 

0.55480E-1 

0.89940E-1 

BIMlI 

— 

0.20770E-1 

0.60834E-1 

0.10024 

giWf* 

— 

0.28109E-1 

0.70127E-1 

0.11052 

0.32991E-2* 

0.35399E-1 

0.78213E-1 

0.12070 

mmWl 

0.51296E-2 

0.42230E-1 

0.88438E-1 

0.13106 

0.69877E-2 

0.49638E-1 

0.97783E-1 

0.14126 

0.125 

■■wiifeMiiaii 

0.68289E-1 

0.12099 

0.16701 

0.25 

0.18327E-1 

0.86652E-1 

0.14465 

0.19269 

0.375 

0.23781E-1 

0.10501 

0.16832 

0.21837 

0.5 

0.29235E-1 

0.12338 

0.19198 

0.24405 

0.625 

0.34690E-1 

0.14174 

0.21563 

0.26835 

0.75 

0.40144E-1 

0.16010 

0.23930 

0.29383 

0.875 

0.17847 

0.26295 

0.31932 

1.0 

0.51053E-1 

0.19683 

0.28468 

0.34959 

IBS 

0.72870E-1 

0.27029 

0.37813 

0.44674 

■w 

0.94087E-1 

0.34375 

0.47158 

0.54870 

3.0 

0.13832 

0.49066 

0.65514 

0.75258 

lEm 

0.63757 

0.84087 

0.95644 

5.0 

0.22559 

0.78448 

1.0265 

1.1603 

*  Gric 

[  independence:  2.2% 

f  Grid  independence:  1.5% 
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Table  9:  Poiseuille  number  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  meniscus  contact  angle  (f3  =  1.0,  0  =  0°). 


Po 

T* 

Tlv 

(j>  —  30° 

'O- 

II 

Oi 

O 

o 

II 

CO 

o 

o 

— 

— 

— 

40.605 

-0.40 

— 

— 

— 

33.803 

-0.35 

— 

— 

40.137 

28.803 

-0.30 

— 

— 

30.719 

25.091 

-0.25 

— 

-0.20 

— 

-0.15 

— 

24.615 

-0.10 

27.918* 

19.546 

mm 

16.384 

0.0 

13.181 

13.939 

13.277 

14.158 

112^ 

10.132 

10.730 

■IlflrIM 

0.25 

5.0257 

7.9851 

8.9752 

0.375 

3.8730 

6.5888 

7.7129 

0.5 

3.1505 

5.6081 

6.7624 

0.625 

2.6551 

4.8816 

wiiwa 

0.75 

2.2944 

0.875 

2.0199 

3.8770 

4.9372 

1.8041 

1.2639 

2.5598 

3.4333 

4.4767 

1 

||l£jQgjQ^fl 

3.6449 

Igjjjpstpi 

2.6575 

4.0 

0.50620 

mm 

IU 

2.0910 

5.0 

0.40829 

0.88203 

1.2646 

1.7236 

*  Gric 

indepenc 

ience:  2.2% 

f  Grid  independence:  1.5% 
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Table  10:  Volumetric  flow  rate  versus  shear  stress  at  the  liquid-vapor  interface  for 
various  values  of  meniscus  contact  angle  (/5  =  1.0,  9  =  0°). 


y* 

Tlv 

o 

O 

II 

~<s>- 

II 

OO 

o 

o 

o 

O 

CO 

II 

-e- 

II 

CD 

O 

o 

— 

— 

— 

0.98398E-1 

— 

— 

— 

0.11832 

e m 

— 

— 

0.52969E-1 

0.13887 

-0.30 

— 

— 

0.69210E-1 

0.15941 

E^i 

— 

— 

0.90857E-1 

0.17988 

-0.20 

— 

0.24531E-1 

0.99625E-1 

0.20048 

— 

0.33199E-1 

0.11484 

0.22104 

gilftl 

0.14160E-2* 

0.41809E-1 

0.12808 

0.24140 

mu 

0.22016E-2 

0.49877E-1 

0.14483 

0.26212” 

0.29991E-2 

0.58626E-1 

0.16013 

0.28252 

0.125 

0.80339E-1 

0.19813 

0.33402 

0.25 

0.78660E-2 

0.10193 

0.23688 

0.38538 

0.375 

0.10270E-1 

0.12353 

0.27565 

0.43674 

0.5 

0.12548E-1 

0.14514 

0.31439 

0.48810 

0.625 

0.14889E-1 

0.16674 

0.35312 

0.53670 

0.75 

0.17230E-1 

0.18834 

0.39189 

0.58766 

0.875 

0.19571E-1 

0.20995 

0.43062 

0.63864 

0.21912E-1 

0.23155 

0.46620 

0.69918 

In 

0.31276E-1 

0.31923 

0.61924 

0.89348 

2.0 

0.40640E-1 

0.40599 

0.77228 

3.0 

0.59368E-1 

0.57951 

■MIS 

4.0 

0.78096E-1 

0.75302 

1.9128 

5.0 

0.96824E-1 

0.92653 

*  Grid  independence:  2.2% 
|  Grid  independence:  1.5% 
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Table  11:  Mean  velocity  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  meniscus  contact  angle  (/3  =  1.0,  6  =  30°). 


v* 

T* 

Tlv 

o 

O 

II 

<f>  =  20° 

It 

O 

o 

II 

05 

O 

o 

-0.375 

— 

— 

— 

0.74804E-1 

-0.3125 

— 

— 

0.41206E-1 

0.93526E-1 

-0.25 

— 

— 

0.56119E-1 

0.11224 

-0.1875 

— 

0.21137E-1 

0.71387E-1 

0.13096 

-0.125 

— 

0.32001E-1 

0.86165E-1 

0.14969 

■iliESEi 

0.57800E-2 

0.42814E-1 

0.10150 

0.16841 

■ 

0.11068E-1 

0.53861E-1 

0.11664 

0.18713 

0.125 

0.21608E-1 

0.74276E-1 

0.14770 

0.22457 

0.25 

0.32361E-1 

0.97689E-1 

0.17876 

0.26202 

0.375 

0.43114E-1 

0.12013 

0.20982 

0.29946 

0.5 

0.53599E-1 

0.14258 

0.23959 

0.33509 

0.625 

0.64249E-1 

0.16503 

0.27028 

0.37225 

0.75 

0.74898E-1 

0.18748 

0.30097 

0.40941 

0.875 

0.85547E-1 

0.20993 

0.33165 

0.44657 

0.96197E-1 

0.23119 

0.36234 

0.48374 

1.5 

0.13879 

0.32003 

0.48508 

0.63238 

0.18139 

0.40887 

0.60782 

0.78103 

0.26658 

0.58655 

0.85008 

0.35177 

0.76423 

1.0943 

■KUM 

0.43697 

0.94191 

1.3386 
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Table  12:  Poiseuille  number  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  meniscus  contact  angle  (/3  =  1.0,  6  =  30°). 


Po 

Tlv 

o 

O 

II 

-e- 

11 

to 

o 

o 

II 

o 

o 

II 

Ob 

o 

o 

-0.375 

— 

— 

— 

-0.3125 

— 

— 

41.349 

-0.25 

— 

-0.1875 

— 

-0.125 

-0.0625 

MSKEm 

0.0 

■HcfkM 

0.125 

MSMM 

0.25 

0.375 

2.9121 

6.4246 

8.1203 

0.5 

2.3424 

5.4127 

yy^ii 

0.625 

1.9541 

4.6764 

6.3039 

0.75 

1.6763 

4.1165 

0.875 

1.4676 

3.6763 

1.3051 

3.3383 

4.7023 

1.5 

0.90460 

2.4116 

0.69217 

1.8876 

0.47096 

1.3158 

0.35691 

1.0098 

0.28732 

0.81939 
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Table  13:  Volumetric  flow  rate  versus  shear  stress  at  the  liquid- vapor  interface  for 
various  values  of  meniscus  contact  angle  (/?  =  1.0,  9  =  30°). 


y* 

Tlv 

^  =  0° 

4>  =  20° 

(j)  =  40° 

4>  =  60° 

-0.375 

— 

— 

— 

0.19279 

-0.3125 

— 

— 

0.81949E-1 

0.24104 

— 

— 

0.14197 

0.28928 

— 

0.28292E-1 

0.11160 

0.33752 

-0.125 

— 

0.42833E-1 

0.17136 

0.38580 

-0.0625 

0.32234E-2 

0.57307E-1 

0.20186 

0.43405 

0.0 

0.59752E-2 

0.72093E-1 

0.23197 

0.48229 

0.125 

0.11665E-1 

0.99419E-1 

0.29374 

0.57879 

0.25 

0.17470E-1 

0.13075 

0.35551 

0.67531 

0.375 

0.23275E-1 

0.16079 

0.41728 

0.77181 

0.5 

0.28936E-1 

0.19084 

0.47651 

0.86364 

0.625 

0.34685E-1 

0.22089 

0.53753 

0.95941 

0.75 

0.40434E-1 

0.25094 

BBliiSlf 

0.875 

0.46183E-1 

0.28099 

0.65959 

1.0 

0.51933E-1 

0.30945 

1.5 

0.75025E-1 

0.42836 

2.0 

0.97926E-1 

0.54727 

1.2088 

3.0 

0.14391 

0.78510 

4.6 

0.18990 

1.0229 

5.0 

I  0.23590 

1.2607 

&xmm 
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Table  14:  Mean  velocity  versus  shear  stress  at  the  liquid-vapor  interface  for  various 
values  of  meniscus  contact  angle  ((3  =  1.0,  6  =  60°). 


V * 

Tlv 

— 

o 

O 

ll 

-e- 

-e- 

II 

H- i 

O 

o 

<f>  =  20° 

II 

CO 

o 

o 

-0.375 

— 

— 

— 

0.75897E-1 

-0.3125 

— 

— 

0.41983E-1 

0.98354E-1 

-0.25 

— 

— 

0.59724E-1 

0.12081 

-0.1875 

— 

0.28621E-1 

0.77466E-1 

0.14326 

-0.125 

0.73396E-2* 

0.41509E-1 

0.94812E-1 

0.16572 

-0.0625 

0.14995E-lt 

0.54181E-1 

0.11263 

0.18818 

0.0 

0.21107 

0.125 

0.37623E-1 

0.92099E-1 

0.16100 

0.25554 

0.25 

0.53077E-1 

0.11829 

0.20175 

0.30046 

0.375 

0.68746E-1 

0.14448 

0.23768 

0.34537 

0.5 

0.84415E-1 

0.17067 

0.27364 

0.39028 

0.625 

0.10008 

0.19648 

0.30959 

0.43296 

0.75 

0.11477 

0.22243 

0.34554 

0.47756 

0.875 

0.13205 

0.24839 

0.38149 

0.52217 

■SB 

0.14776 

0.27435 

0.41744 

0.56677 

1.5 

■MM 

0.37816 

0.56125 

0.74520 

0.26815 

0.48199 

0.70506 

0.92361 

0.39098 

0.68963 

0.51380 

0.89727 

0.63663 

1.1049 

*  Grid  independence:  2.6% 
f  Grid  independence:  1.1% 
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Table  15:  Poiseuille  number  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  meniscus  contact  angle  (/?  =  1.0,  6  —  60°). 


Po 

rlv 

II 

O 

o 

o 

O 

r-i 

II 

-e- 

II 

to 

o 

o 

o 

O 

CO 

II 

-0.375 

— 

— 

— 

40.780 

-0.3125 

— 

— 

43.296 

31.469 

-0.25 

— 

— 

30.434 

25.619 

-0.1875 

— 

30.202 

23.464 

21.604 

-0.125 

31.963* 

20.699 

19.171 

18.676 

-0.0625 

15.644t 

15.858 

16.138 

16.447 

0.0 

10.361 

12.833 

13.932 

14.663 

0.125 

6.2354 

9.3294 

10.943 

12.111 

0.25 

4.4199 

7.2637 

9.0097 

10.301 

0.375 

3.4125 

5.9470 

7.6474 

8.9617 

0.5 

2.7791 

5.0344 

6.6426 

7.9304 

0.625 

2.3440 

4.3730 

5.8712 

7.1487 

0.75 

2.0440 

3.8627 

5.2604 

6.4810 

0.875 

1.8017 

3.4591 

4.7646 

5.9274 

1.0 

1.6108 

3.1318 

4.3543 

5.4609 

1.5 

1.1347 

2.2720 

3.2386 

4.1534 

2.0 

0.87487 

1.7826 

2.5780 

3.3526 

3.0 

0.60003 

1.2459 

1.8311 

2.4172 

4.0 

0.45659 

0.95760 

1.4197 

1.8904 

5.0 

0.36850 

0.77764 

1.1593 

1.5521 

*  Grid  independence:  2.6% 
f  Grid  independence:  1.1% 
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Table  16:  Volumetric  flow  rate  versus  shear  stress  at  the  liquid- vapor  interface  for 
various  values  of  meniscus  contact  angle  (/?  =  1.0,  0  =  60°). 


V* 

II 

o 

o 

<t>  =  10° 

II 

to 

o 

o 

<p  =  30° 

-0.375 

— 

— 

— 

0.28325 

-0.3125 

— 

— 

0.12007 

0.36706 

-0.25 

— 

— 

0.17081 

0.45086 

-0.1875 

— 

0.56279E-1 

0.22155 

0.53465 

-0.125 

0.75407E-2* 

0.81621E-1 

0.27116 

0.61847 

-0.0625 

0.15405E-1T 

0.10651 

0.32212 

0.70229 

0.0 

0.23264E-1 

0.13164 

0.37311 

0.78772 

0.125 

0.38656E-1 

0.18109 

0.46046 

0.95368 

0.25 

0.54535E-1 

0.23260 

0.57700 

1.1213 

0.375 

0.70634E-1 

0.28409 

0.67976 

1.2889 

0.5 

0.86734E-1 

0.33559 

0.78261 

1.4565 

0.625 

0.10282 

0.38635 

0.88543 

1.6158 

0.75 

0.11792 

0.43739 

0.98825 

1.7822 

|  0.875 

0.13567 

0.48843 

1.0910 

1.9487 

1.0 

0.15181 

0.53947 

1.1940 

2.1152 

1.5 

0.21242 

0.74359 

1.6051 

2.7811 

2.0 

0.27551 

0.94776 

2.0164 

3.4469 

3.0 

0.40172 

1.3560 

2.8390 

4.7785 

4.0 

0.52791 

1.7643 

3.6613 

6.1101 

5.0 

0.65412 

2.1726 

4.4842 

7.4420 

*  Grid  independence:  2.6% 
f  Grid  independence:  1.1% 
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Table  17:  Mean  velocity  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  groove  aspect  ratio  ( <j>  =  30°,  0  =  0°). 


V* 

Tlv 

0  =  0.5 

0  =  1.0 

0  =  1.5 

-0.375 

0.14752E-1* 

— 

— 

-0.3125 

0.19853E-1** 

— 

— 

-0.25 

0.25245E-1 

— 

— 

-0.1875 

0.30603E-1 

0.22287E-ln 

— 

-0.125 

0.35952E-1 

0.31469E-1** 

0.12684E-1* 

-0.0625 

0.41301E-1 

0.40743E-1 

0.19402E-1 

0.0 

0.45883E-1 

0.49638E-1 

0.26119E-1 

0.125 

0.56930E-1 

0.68289E-1 

0.39552E-1 

0.25 

0.67443E-1 

0.86652E-1 

0.52986E-1 

0.375 

0.77955E-1 

0.10501 

0.66420E-1 

0.5 

0.88468E-1 

0.12338 

0.79853E-1 

0.625 

0.98981E-1 

0.14174 

0.93287E-1 

0.75  1 

0.10949 

0.16010 

0.10672 

0.875 

0.12047 

0.17847 

0.12015 

0.13102 

0.19683 

0.13358 

1.5 

0.17326 

0.27029 

0.18732 

0.24105 

o.sooosf 

0.49066 

0.34852 

0.38444** 

0.63757 

0.45599 

5.0 

0.46891** 

0.78448 

0.56346 

*  Grid  independence:  2.2% 
**  Grid  independence:  1.2% 
f  Grid  independence:  1.1% 
ft  Grid  independence:  2.0% 
\  Grid  independence:  2.3% 
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Table  18:  Poiseuille  number  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  groove  aspect  ratio  (<fi  =  30°,  9  —  0°). 


Po 

Tlv 

f3  =  0.5 

o 

r“H 

II 

II 

Cn 

-0.375 

38.108* 

— 

— 

-0.3125 

28.317** 

— 

— 

-0.25 

22.260 

— 

— 

-0.1875 

18.370 

31.045** 

— 

-0.125 

15.636 

21.987** 

33.179* 

-0.0625 

13.611 

16.982 

21.692 

0.6 

11.976 

13.939 

16.114 

0.125 

9.8457 

10.132 

10.641 

0.25 

8.3110 

7.9851 

7.9434 

0.375 

7.1902 

6.5888 

6.3368 

0.5 

6.3358 

5.6081 

5.2707 

0.625 

5.6629 

4.8816 

4.5117 

1  0.75 

5.1192 

4.3217 

3.9438 

0.875 

4.6665 

3.8770 

3.5029 

1.0 

4.2905 

3.5153 

3.1506 

1.5 

3.2446 

2.5598 

2.2468 

2.0 

2.6093* 

2.0128 

1.7460 

3.0 

1.8736* 

1.4102 

1.2076 

4.0 

1.4623** 

1.0852 

0.92301 

5.0 

1.1989** 

0.88203 

0.74697 

*  Grid  independence:  2.2% 
**  Grid  independence:  1.2% 
f  Grid  independence:  1.1% 
ff  Grid  independence:  2.0% 
X  Grid  independence:  2.3% 
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Table  19:  Volumetric  flow  rate  versus  shear  stress  at  the  liquid- vapor  interface  for 
various  values  of  groove  aspect  ratio  ((f)  =  30°,  0  =  0°). 


y* 

T* 

Tlv 

f3  =  0.5 

13  =  1.0 

(3  =  1.5 

-0.375 

0.11731E-1* 

— 

— 

-0.3125 

0.15788E-1** 

— 

— 

-0.25 

0.20084E-1 

— 

— 

-0.1875 

0.24337E-1 

0.26218E-1TT 

— 

-0.125 

0.28591E-1 

0.37020E-1** 

0.14548E-1* 

-0.0625 

0.32845E-1 

0.47929E-1 

0.22456E-1 

0.0 

0.36076E-1 

0.58626E-1 

0.29954E-1 

0.125 

0.45208E-1 

0.80339E-1 

0.45361E-1 

0.25 

0.53556E-1 

0.10193 

0.60767E-1 

0.375 

0.61904E-1 

0.12353 

0.76174E-1 

0.5 

0.70252E-1 

0.14514 

0.91580E-1 

0.625 

0.78600E-1 

0.16674 

0.10698 

0.75 

0.86948E-1 

0.18834 

0.12239 

i  0.875 

0.95806E-1 

0.20995 

0.13780 

1.0 

0.10420 

0.23155 

0.15320 

1.5 

0.13779 

0.31923 

0.21483 

2.0 

0.17134* 

0.40599 

0.27645 

3.0 

0.23862+ 

0.57951 

0.39971 

4.0 

0.30573** 

0.75302 

0.52296 

5.0 

0.37291** 

0.92653 

0.64621 

*  Grid  independence:  2.2% 
**  Grid  independence:  1.2% 
f  Grid  independence:  1.1% 
ft  Grid  independence:  2.0% 
|  Grid  independence:  2.3% 


101 


Table  20:  Mean  velocity  versus  shear  stress  at  the  liquid-vapor  interface  for  various 
values  of  groove  aspect  ratio  ( <fr  =  30°,  6  =  30°). 


V* 

Tlv 

O 

o 

II 

(3  =  0.5 

(3  =  1.0 

0  =  1.5 

-0.25 

— 

0.33511E-1 

0.31252E-1 

— 

-0.1875 

0.10930E-1 

0.44442E-1 

0.44545E-1 

0.30859E-1* 

-0.125 

0.15939E-1 

0.55373E-1 

0.57663E-1 

0.43708E-lf 

-0.0625 

0.21027E-1 

0.66304E-1 

0.70993E-1 

0.56557E-1 

0.0 

0.26234E-1 

0.77458E-1 

0.83936E-1 

0.68965E-1 

0.125 

0.36293E-1 

0.99459E-1 

0.11017 

0.94766E-1 

0.25 

0.46235E-1 

0.12154 

0.13723 

0.12056 

0.375 

0.56325E-1 

0.14362 

0.16429 

0.14636 

0.5 

0.66416E-1 

0.16468 

0.19135 

0.17216 

0.625 

0.76507E-1 

0.18654 

0.21841 

0.19796 

0.75 

0.86598E-1 

0.20840 

0.24547 

0.22376 

0.875 

0.96689E-1 

0.23026 

0.27253 

0.24956 

0.10678 

0.25212 

0.29959 

0.27537 

1.5 

0.14714 

0.33957 

0.40545 

0.37857 

0.18750 

0.42702 

0.51262 

0.48177 

0.26823 

0.60191 

0.72697 

0.68818 

0.34896 

0.77680 

0.94131 

0.89458 

0.42968 

0.95170 

1.1556 

1.1009 

*  Grid  independence:  1.9% 
f  Grid  independence:  1.3% 
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Table  21:  Poiseuille  number  versus  shear  stress  at  the  liquid- vapor  interface  for  various 
values  of  groove  aspect  ratio  (<f)  =  30°,  6  —  30°). 


Po 

Tlv 

O 

o 

II 

/?  =  0.5 

13  =  1.0 

(3  =  1.5 

-0.25 

— 

29.167 

38.710 

— 

-0.1875 

28.606 

21.993 

27.158 

37.289* 

-0.125 

19.617 

17.651 

20.980 

26.3271' 

-0.0625 

14.869 

14.741 

17.040 

20.346 

0.0 

11.918 

12.618 

14.413 

16.686 

0.125 

8.6153 

9.8276 

10.980 

12.143 

0.25 

6.7628 

8.0422 

8.8153 

9.5446 

0.375 

5.5513 

6.8057 

7.3635 

7.8621 

0.5 

4.7078 

5.9353 

6.3222 

6.6839 

0.625 

4.0869 

5.2398 

5.5390 

5.8126 

0.75 

3.6107 

4.6901 

4.9284 

5.1426 

0.875 

3.2338 

4.2448 

4.4391 

4.6109 

1.0 

2.9282 

3.8767 

4.0381 

4.1789 

1.5 

2.1250 

2.8784 

2.9838 

3.0397 

2.0 

1.6675 

2.2889 

2.3600 

2.3885 

3.0 

1.1657 

1.6239 

1.6641 

1.6721 

4.0 

0.89603 

1.2582 

1.2852 

1.2863 

5.0 

0.72769 

1.0270 

1.0468 

1.0452 

*  Grid  independence:  1.9% 
f  Grid  independence:  1.3% 
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Table  22:  Volumetric  flow  rate  versus  shear  stress  at  the  liquid-vapor  interface  for 
various  values  of  groove  aspect  ratio  {(j>  =  30°,  6  =  30°). 


V* 

■HI 

o 

o 

II 

/3  =  0.5 

(3  =  1.0 

(3  =  1.5 

— 

0.38764E-1 

0.52371E-1 

— 

mmm 

0.49904E-2 

0.51410E-1 

0.74649E-1 

0.62142E-1* 

-0.125 

0.72773E-2 

0.64055E-1 

0.96632E-1 

o.ssoieE-i1 

-0.0625 

0.96006E-2 

0.76699E-1 

0.11897 

0.11389 

0.0 

0.11977E-1 

0.89602E-1 

0.14066 

0.13887 

0.125 

0.16570E-1 

0.11505 

0.18464 

0.19083 

0.25 

0.21109E-1 

0.14059 

0.22998 

0.24278 

0.375 

0.25716E-1 

0.16613 

0.27532 

0.29472 

0.5 

0.30323E-1 

0.19050 

0.32066 

0.34668 

0.625 

0.34931E-1 

0.21578 

0.36602 

0.39846 

0.75 

0.39538E-1 

0.24106 

0.41136 

0.45060 

0.875 

0.44145E-1 

0.26636 

0.45670 

0.50245 

0.48752E-1 

0.29164 

0.50206 

0.55450 

1.5 

0.67181E-1 

0.39280 

0.67946 

0.76232 

0.85609E-1 

0.49396 

0.85906 

0.97014 

0.12246 

0.69628 

1.2182 

MiliiUkVJM 

0.89859 

1.5774 

5.0 

0.19618 

1.1009 

1.9366 

*  Grid  independence:  1.9% 
f  Grid  independence:  1.3% 
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Table  23:  Mean  velocity  versus  shear  stress  at  the  liquid- vapor  interface  for 
values  of  groove  aspect  ratio  (<j>  =  30°,  6  =  60°). 


V * 

T* 

rlv 

o 

o 

II 

/?  =  0.5 

ft  — 1.0 

p  =  1.5 

-0.375 

— 

0.53334E-1 

0.75897E-1 

0.90127E-1 

-0.3125 

0.33653E-1 

0.72695E-1 

0.98354E-1 

0.11443 

-0.25 

0.48754E-1 

0.92056E-1 

0.12081 

0.13874 

-0.1875 

0.63856E-1 

0.11141 

0.14326 

0.16306 

-0.125 

0.78957E-1 

0.13077 

0.16572 

0.18737 

-0.0625 

0.94058E-1 

0.15013 

0.18818 

0.21168 

0.0 

0.10915 

0.17030 

.0.21107 

0.23708 

0.125 

0.13936 

0.20824 

0.25554 

0.28461 

0.25 

0.16956 

0.24694 

0.30046 

0.33324 

0.375 

0.19976 

0.28566 

0.34537 

0.38185 

0.5 

0.22996 

0.32438 

0.39028 

0.43047 

0.625 

0.26017 

0.36311 

0.43296 

0.47910 

0.75 

0.29037 

0.40183 

0.47756 

0.52772 

0.875 

0.32057 

0.44055 

0.52217 

0.57634 

1.0 

0.35080 

0.47927 

0.56677 

0.62496 

0.47161 

0.63416 

0.74520 

0.81945 

0.58941 

0.78905 

0.92361 

1.0088" 

0.82971 

1.0988 

■11HH 

1.3955 

IBESSifll 

2.1688 

various 
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Table  24:  Poiseuille  number  versus  shear  stress  at  the  liquid-vapor  interface  for  various 
values  of  groove  aspect  ratio  (<£  =  30°,  6  =  60°). 


Po 

Tlv 

(3  =  0.0 

(3  =  0.5 

Xb 

II 

(— 1 

o 

(3=  1.5 

-0.375 

— 

44.783 

40.780 

-0.3125 

44.571 

32.856 

31.469 

-0.25 

30.766 

25.946 

25.619 

-0.1875 

23.490 

21.437 

21.604 

-0.125 

■HEEHi 

18.676 

-0.0625 

16.447 

HI 

0.0 

13.741 

14.025 

14.663 

0.125 

KH 

0.25 

mb 

BOSS 

mmm 

0.375 

8.9617 

0.5 

0.625 

0.75 

ppf.sgi 

MESS 

0.875 

5.4216 

ual 

4.2758 

4.9835 

5.4609 

1.5 

glBfeSM 

2.5448 

KBiEHSB 

|  | 

■*:«)rWi 

wtimsm 

1.1447 

1.3899 

1.5521 
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Table  25:  Volumetric  flow  rate  versus  shear  stress  at  the  liquid- vapor  interface  for 
various  values  of  groove  aspect  ratio  (<f>  =  30°,  9  =  60°). 


y* 

T* 

Tlv 

II 

o 

o 

£  =  0.5 

O 

i—l 

II 

£  =  1.5 

-0.375 

— 

0.14571 

0.28325 

0.42648 

-0.3125 

0.58288E-1 

0.19860 

0.36706 

0.54152 

-0.25 

0.84444E-1 

0.25150 

0.45086 

0.65656 

-0.1875 

0.11060 

0.30439 

0.53465 

0.77160 

-0.125 

0.13675 

0.35729 

0.61847 

0.88664 

-0.0625 

0.16291 

0.41018 

0.70229 

1.0016 

0.0 

0.18907 

0.46526 

0.78772 

1.1219 

0.125 

0.24138 

0.56892 

0.95368 

1.3468 

0.25 

0.29368 

0.67466 

1. 1213m 

1.5769 

0.375 

0.34600 

0.78045 

1.2889 

1.8069 

0.5 

0.39830 

0.88624] 

1.4565 

2.0370 

0.625 

0.45062 

0.99203 

1.6158 

2.2670 

0.75 

0.50294 

1.0978 

1.7822 

2.4972 

0.875 

0.55524 

1.2036 

1.9487 

2.7272 

1.0 

0.60760 

1.3094 

2.1152 

2.9573 

1.5 

0.81686 

1.7325 

2.7811 

3.8776 

2.0 

1.0208 

2.1557 

3.4469 

4.7737 

3.0 

1.4370 

3.0020 

4.7785 

6.6035 

4.0 

1.8532 

3.8483 

6.1101 

8.4334 

5.0 

2.2695 

4.6947 

7.4420 

10.263 
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Table  26:  Liquid-vapor  shear  stress  when  the  mean  velocity  is  zero  versus  meniscus 
contact  angle  for  f3  =  1.0. 


~Tlv.  0 

Analytical 

Numerical 

<t> 

II 

O 

o 

o 

O 

CO 

II 

0  =  60° 

II 

Co 

o 

o 

o 

O 

CO 

II 

0° 

0.13662 

0.14152 

0.17957 

0.16743 

0.12910 

0.18826 

1 

o 

o 

— 

— 

0.35260 

— 

— 

0.32225 

20° 

— 

0.39066 

0.52077 

— 

0.30236 

0.45252 

o 

O 

CO 

0.48838 

— 

— 

0.59013 

o 

O 

— 

0.61770 

— 

— 

0.48096 

— 

o 

o 

0.78192 

0.81699  ■ 

— 

0.53204 

0.62874 

— 

90° 

1.0000 

— 

— 

0.692981 

— 

— 

Table  27:  Liquid-vapor  shear  stress  when  the  mean  velocity  is  zero  versus  groove 
aspect  ratio  for  <j>  =  30°. 


~Tlv.O 

Analytical 

Numerical 

0 

II 

O 

o 

o 

O 

CO 

II 

6  =  60° 

II 

O 

o 

II 

CO 

o 

o 

o 

O 

CO 

II 

OS 

0 

— 

0.37758 

0.50000 

— 

0.32342 

0.45198 

0.5 

0.65768 

0.51267 

0.61200 

0.54363 

0.44253 

0.54762 

1 

0.48838 

0.50727 

0.39337 

0.59013 

1.5 

0.31907 

0.46283 

0.73205 

0.24279 

0.33430 

0.61283 
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Table  28:  Mean  velocity,  Poiseuille  number  and  volumetric  flow  rate  versus  groove 
half-angle  {<j>  =  30.0°,  r,;  =  5.0). 


0 

e 

V * 

Po 

V* 

0.0 

15° 

0.15039 

0.63916 

0.34134E-1 

0.0 

CO 

o 

o 

0.42968 

0.72769 

0.19618 

0.0 

45° 

0.77729 

0.87321 

0.64037 

0.0 

60° 

1.3103 

1.1447 

2.2695 

0.5 

0° 

0.46840* 

1.2002* 

0.37250* 

0.5 

15° 

0.72524 

1.0147 

0.88294 

0.5 

o 

O 

CO 

0.95170 

1.0270 

1.1009 

0.5 

45° 

1.2307 

1.1405 

1.9737 

0.5 

60° 

1.7183 

1.3899 

4.6947 

1.0 

0° 

0.78448 

0.88203 

0.92653 

1.0 

15° 

0.95403 

0.92271 

1.2882 

1.0 

o 

O 

CO 

1.1556 

1.0468 

1.9366 

1.0 

45° 

1.4570 

1.2408 

3.3446 

1.0 

60° 

1.9941 

1.5521 

7.4420 

1.5 

0° 

0.56346 

0.74697 

0.64621 

1.5 

15° 

0.83014 

0.82527 

1.2318 

1.5 

30° 

1.1009 

1.0452 

2.2170 

1.5 

45° 

1.15058 

1.3144 

4.3656 

1.5 

60° 

2.1688 

1.6856 

10.263 

*  Grid  independence:  1.9% 
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Table  29:  Mean  velocity,  Poiseuille  number  and  volumetric  flow  rate  versus  groove 
half-angle  (/3  =  1.0,  r;*  =  5.0). 


4> 

e 

v* 

Po 

y* 

0° 

0° 

0.22559 

0.40829 

0.96824E-1 

0° 

15° 

0.35819 

0.26472 

0.15874 

0° 

GO 

O 

o 

0.43697 

0.28732 

0.23590 

0° 

45° 

0.51700 

0.34103 

0.37059 

0°  1 

60° 

0.63663 

0.36850 

0.65412 

10° 

0° 

0.41271* 

0.65741* 

0.30402* 

10° 

15° 

0.56551+ 

0.54255+ 

0.45080t 

10° 

o 

O 

CO 

0.69597 

0.57652 

0.67168 

10° 

45° 

0.85547 

0.66278 

1.0996 

10° 

60° 

1.1049 

0.77764 

2.1726 

20° 

0° 

0.63489 

0.75699 

0.62239 

o 

o 

CN 

15° 

0.79550 

0.72604 

0.87007 

20° 

o 

O 

CO 

0.94191 

0.81939 

1.2607 

to 

o 

o 

45° 

1.1689 

0.95753 

2.1112 

20° 

60° 

1.5679 

1.1593 

4.4842 

o 

O 

CO 

0° 

0.78448 

0.88203 

0.92653 

o 

O 

CO 

15° 

0.95512 

0.92166 

1.2896 

o 

O 

CO 

o 

O 

CO 

1.1556 

1.0468 

1.9366 

GO 

O 

o 

45° 

1.4570 

1.2408 

3.3446 

o 

O 

CO 

60° 

1.9941 

1.5521 

7.4420 

o 

O 

0° 

0.91623 

0.99751 

1.2387 

o 

O 

15° 

1.0914 

1.1027 

1.7231 

o 

O 

o 

O 

CO 

1.3386 

1.2727 

2.6621 

o 

O 

45° 

1.7251 

1.5229 

4.7736 

Cn 

O 

o 

0° 

0.98256 

1.1480 

1.4758 

50° 

15° 

1.2066 

1.2793 

2.1575 

50° 

o 

O 

CO 

1.5161 

1.4856 

3.4670 

60° 

0° 

1.0265 

1.2646 

1.6810 

60° 

15° 

1.3069 

1.4552 

2.5940 

cn> 

o 

o 

o 

O 

CO 

1.6723 

1.7105 

4.3101 

*  Grid  independence:  1.5% 
f  Grid  independence:  1.7% 
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Table  30:  Geometric  values  for  the  parametric  analysis  to  determine  the  volumetric 
flow  rate  versus  groove  fill  ratio  ( h  =  3.831  x  10-4  m,  w  =  3.445  x  10-4  m,  4> o  =  7.0°, 
Ag  =  1.703  x  1(T7  m2). 


p 

e 

ht  (m) 

Wi  (m) 

R  (m) 

Ai  (m2) 

Al/Ag 

0.0 

52.31° 

7.0° 

1.653  x  10“° 

4.279x10-° 

4.192x10"° 

3.674xl0_lu 

0.00216 

0.0 

52.31° 

7.0° 

3.401  x  10"° 

8.804X10-5 

8.625x10-° 

1.555xl0-y 

0.00913 

0.0 

52.31° 

7.0° 

5.140  x  lO"5 

1.331  xlO"4 

1.303xl0-4 

3.552xl0"y 

0.02086 

0.0 

52.31° 

7.0° 

6.879  x  10“° 

1.781  xlO"4 

1.744xlO-4 

6.362xl0-9 

0.03737 

1.049 

14.62° 

7.0° 

1.642  x  10~4 

4.302  xlO"4 

2.314xl0-4 

1.805x10-° 

0.10603 

0.8991 

14.62° 

7.0° 

1.916  x  10“4 

4.445  xlO"4 

2.390xl0"4 

2.696x10-° 

0.15833 

0.7868 

14.62° 

7.0° 

2.189  x  lO"4 

4.587  xlO-4 

2.467xlO-4 

3.612x10"° 

0.21214 

0.6993 

14.62° 

7.0° 

2.463  x  10"4 

4.730xl0"4 

2.544xlO-4 

4.561x10-° 

0.26785 

0.6294 

14.62° 

7.0° 

2.737  x  10“4 

4.873  xlO"4 

2.621  xlO"4 

5.538x10-® 

0.32528 

0.5722 

14.62° 

7.0° 

3.010  x  10"4 

5.015xl0-4 

2.697xl0-4 

6.541x10"° 

0.38419 

0.5245 

14.62° 

7.0° 

3.284  x  10~4 

5.158xl0“4 

2.774xlO"4 

7.577x10-° 

0.44503 

0.4841 

14.62° 

7.0° 

3.558  x  10~4 

5.301  xlO"4 

2.851X10"4 

8.642x10“° 

0.50757 

0.4496 

14.62° 

7.0° 

3.831  x  lO"4 

5.444xl0-4 

2.928xl0-4 

9.732x10“° 

0.57159 

0.4496 

14.62° 

15.0° 

3.831  x  10“4 

5.444  xlO-4 

3.131  xlO"4 

1.091  xlO"7 

0.64063 

0.4496 

14.62° 

25.0° 

3.831  x  10"4 

5.444  xlO"4 

3.534xlO"4 

1.218xlO_Y 

0.71541 

0.4496 

14.62° 

35.0° 

3.831  x  10~4 

5.444xlO-4 

4.201  xlO"4 

1.330  xlO~Y 

0.78100 

0.4496 

14.62° 

45.0° 

3.831  x  10"4 

5.444xl0-4 

5.382xl0-4 

1.431xl0"7 

0.84019 

0.4496 

14.62° 

55.0° 

3.831  x  10"4 

5.444  xlO"4 

7.816xlO-4 

1.524xlO-Y 

0.89505 

0.4496 

14.62° 

65.0° 

3.831  x  10~4 

5.444xl0-4 

1.511xlO-3 

1.613xl0-v 

0.94722 

0.4496 

14.62° 

75.0° 

3.831  x  10“4 

5.444xlO-4 

4.104X10-3 

1.699  xlO"7"" 

0.99808 

Ill 


Table  31:  Mean  velocity,  volumetric  flow  rate  and  dimensional  flow  parameters  versus 
groove  fill  ratio  ( h  =  3.831  x  10-4  m,  w  =  3.445  x  10-4  m,  0o  =  7.0°,  t*v  =  0.0). 


A, /A, 

^0* 

V* 

Fi  (m2) 

F2  (m4) 

0.00216 

0.35025X10'1 

0.23545  xlO"1 

rnmria 

0.00913 

0.35025xl0_1 

0.23545  xlO"1 

0.02086 

0.35025xl0_1 

0.23545  xlO'1 

Hsmaiai 

m&MivmiimM 

0.03737 

0.35025xl0_1 

0.23545xl0-1 

1.0545xl0“18 

■nwraiKl 

0.15833 

0.20927xl0_1 

0.15367xl0-1 

witrmiia 

0.21214 

0.27194xl0_1 

0.20498  xlO"1 

0.26785 

0.32454  xlO-1 

0.24398xl0_1 

rrmoraai 

0.32528 

0.36072  xlO'1 

0.26669X10"1 

0.38419 

0.38622xl0_1 

0.27884xl0_1 

g 

0.44503 

0.39823xl0_1 

0.27979X10"1 

0.50757 

0.40251  xlO"1 

0.27477xl0-1 

5.0955  xlO"9 

4.4035  xlO" 16 

0.57159 

0.40053X10-1 

msmEM 

K2stts£Si!fiil 

5.7208  xlO'16 

0.64063 

0.46394  xlO"1 

0.34479  xlO"1 

MraaSEWiijjg 

0.71541 

0.53536xl0_1 

0.44431  xlO-1 

0.78100 

0.59941X10"1 

0.54307xl0_1 

1.1698x10“ 15 

gSggggggigjg 

0.89505 

0.70861  xlO"1 

0.73575  xlO"1 

[Em 

0.94722 

0.75803xl0_1 

0.83295  xlO"1 

0.99808 

0.80163xl0_1 

0.92993  xlO-1 

2.0031  xlO"15 
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Table  32:  Dimensionless  velocity,  volumetric  flow  rate  and  Poiseuille  number  versus 
groove  fill  ratio  ( h  =  3.831  x  10-4  m,  w  =  3.445  x  10-4  m,  4>q  =  7.0°,  7}*  a  —  5.0). 


0 

9 

<t> 

A  l/Ag 

K 

V*"~ 

Po 

1.049 

14.62° 

7.0° 

0.10603 

0.44951 

0.45997 

0.30098 

0.8991 

14.62° 

7.0° 

0.15833 

0.55128 

0.52375 

0.40480 

0.7868 

14.62° 

7.0° 

0.21214 

0.61118* 

0.46070* 

0.56117* 

0.6993 

14.62° 

7.0° 

0.26785 

0.65735** 

0.49417** 

0.57269** 

0.6294 

14.62° 

7.0° 

0.32528 

0.66535 

0.59419 

0.49190 

0.5722 

14.62° 

7.0° 

0.38419 

0.66313 

0.60976 

0.47876 

0.5245 

14.62° 

7.0°  1 

0.44503 

0.64877** 

0.45581** 

0.62692** 

0.4841 

14.62° 

7.0° 

0.50757 

0.62513 

0.64739 

0.42674 

0.4496 

14.62° 

7.0° 

0.57159 

0.59685 

0.66990 

0.39574 

0.4496 

14.62° 

15.0° 

0.64063 

0.62644 

0.80169 

0.46555 

0.4496 

14.62° 

25.0° 

0.71541 

0.65255 

0.95976 

0.54157 

0.4496 

14.62° 

35.0° 

0.78100 

0.66745 

1.1183 

0.60472 

0.4496 

14.62° 

45.0° 

0.84019 

0.68388 

1.2631 

0.66656 

0.4496 

14.62° 

55.0° 

0.89505 

0.69685 

1.4068 

0.72354 

0.4496 

14.62° 

65.0° 

0.94722 

0.71080 

1.5446 

0.78105 

0.4496 

14.62° 

75.0° 

0.99808 

0.72674 

1.6838 

0.84305 

*  Grid  independence:  1.5% 
**  Grid  independence:  1.2% 


Table  33:  Maximum  heat  transport  predicted  by  the  closed-form  solution  versus 
groove  fill  ratio  (Straight  axial  grooves,  no  body  forces,  Tsat  =  40°C,  Ethanol) 


Q«p  (W) 

A, /A, 

Semi-Analytical 

Two-Point  Numerical 

No  Shear 

0.10598 

0.13573 

0.13570 

0.13590 

0.15830 

0.39121 

0.39097 

0.39217 

0.21209 

0.88732 

0.88654 

0.89115 

0.26782 

1.6887 

1.6867 

0.32519 

2.8028 

0.38408 

4.2781 

4.2732 

0.44492 

6.0615 

6.0553 

0.50745 

8.1727 

8.1669 

0.57146 

10.575 

0.64063 

13.703 

13.676 

0.71520 

17.609 

17.541 

0.78097 

21.483 

21.369 

be m 

0.84028 

25.247 

25.072 

0.89489 

28.992 

28.753 

0.94715 

32.759 

32.449 

0.99765 

36.416 

36.023 
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Table  34:  Maximum  heat  transport  predicted  by  the  closed-form  solution  versus 
groove  fill  ratio  (Straight  axial  grooves,  no  body  forces,  Tsa t  =  40°C,  Water) 


CU(W) 

MBEM 

Semi- Analytical 

Two-Point  Numerical 

No  Shear 

0.10598 

1.7825 

1.7787 

1.8044 

0.15830 

5.0811 

5.0508 

5.2059 

0.21209 

11.338 

11.243 

0.26782 

0.32519 

34.292 

33.915 

gaii 

0.38408 

50.924 

50.408 

0.44492 

70.242 

69.617 

0.50745 

92.136 

91.584 

0.57146 

116.02 

115.79 

0.64063 

147.75 

145.45 

0.71520 

186.60 

181.02 

0.78097 

224.33 

215.29 

0.84028 

260.28 

246.91 

0.89489 

295.11 

277.51 

0.94715 

329.09 

307.01 

0.99765 

361.03 

333.96 
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Table  35:  Maximum  heat  transport  versus  groove  fill  ratio  for  several  working  tem¬ 
peratures  (Ethanol,  |ar|  =  0.0-g). 


Qm,  (W) 

AJA, 

Tsat  =  40°C 

Tsat  =  60°  C 

P» 

II 

oo 

CD 

o 

o 

Tsat  =  100°C 

0.10603 

0.0 

0.0 

0.12 

0.15833 

6.20 

0.35 

1.13 

0.21214 

0.68 

2.82 

0.26785 

1.47 

3.26 

5.55 

0.32528 

2.60 

3.59 

5.53 

0.38419 

IKBHi 

8.58 

0.44503 

5.86 

12.25 

NEliSNI 

0.50757 

7.89 

10.98 

16.66 

0.57159 

0.64063 

■n 

18.63 

28.12 

47.21 

0.71541 

17.43 

24.03 

36.34 

0.78100 

29.40 

44.42 

0.84019 

25.06 

34.59 

52.29 

87.63 

0.89505 

28.85 

39.79 

0.94722 

32.65 

45.02 

0.99808 

36.32 

50.13 

75.85 

127.22 
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Table  36:  Maximum  heat  transport  versus  groove  fill  ratio  for  several  working  tem¬ 
peratures  (Ethanol,  \ar\  —  10.0-g). 


Qcap  (W) 

Al/Ag 

Tsat  =  40°  C 

Tsat  =  60°C 

Tsat  =  80°  C 

Tsat  =  100°C 

0.10603 

0.34 

0.59 

0.99 

1.87 

0.15833 

1.29 

1.92 

3.11 

5.63 

3.11 

4.53 

7.24 

12.95 

0.26785 

6.07 

8.73 

0.32528 

10.20 

14.60 

41.44 

0.38419 

22.44 

35.66 

63.40 

0.44503 

50.72 

90.16 

0.50757 

30.00 

68.63 

121.94 

0.57159 

38.97 

56.05 

89.13 

158.41 

0.64063 

50.55 

115.67 

0.71541 

93.73 

149.02 

114.49 

0.84019 

93.52 

134.75 

| 

0.89505 

107.55 

155.01 

246.58 

438.41 

0.94722 

121.71 

279.05 

496.46 

0.99808 

135.45 

195.30 

310.81 

552.75 
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Table  37:  Maximum  heat  transport  versus  radial  acceleration  ( Ai/Ag  =  0.5). 


<3«p  (W) 

Castle  et  al.  (2000) 

Present 

£1  (g) 

No  Dryout 

Partial  Dryout 

Analytical 

Analytical 

0.01 

— 

— 

19 

7.40 

2.0 

— 

— 

22 

9.56 

4.0 

— 

— 

31 

13.78 

6.0 

42 

18.32 

8.0 

9.5 

13.0 

58 

22.95 

10.0 

11.0 

13.5 

70 

27.65 

Table  38:  Maximum  heat  transport  versus  radial  acceleration  ( Ai/Ag  =  1.0). 


<3« p  (W) 

Castle  et  al.  (2000) 

Present 

1*1  (g) 

No  Dryout 

Partial  Dryout 

Analytical 

Analytical 

0.01 

17.3 

25.4 

20 

34.97 

2.0 

17.3 

25.4 

24 

44.64 

4.0 

32.7 

40.0 

39 

70.29 

6.0 

59.1 

64.5 

52 

98.07 

8.0 

88.2 

96.4 

75 

138.04 

10.0 

106.4 

110.9 

92 

172.67 

118 


Table  39:  Mean  velocity,  Poiseuille  number  and  volumetric  flow  rate  versus  shear 
stress  at  the  liquid-vapor  interface  for  various  values  of  meniscus  contact  angle  (/?  = 
0.5,  wf/2  =  0.25,  P*  =  1.15245). 


<t> 

Tb>, 

Po 

V* 

0° 

-0.1 

0.18878  x  10-2- 

27.088 

0.17394  x  10-° 

0° 

-0.075 

0.26896  x  10~* 

19.013 

0.24782  x  10“° 

0° 

-0.05 

0.35104  x  io_* 

14.567 

0.32344  x  10“° 

0° 

-0.025 

0.43465  x  10"2 

11.765 

0.40048  x  10“° 

0° 

0.0 

0.51353  x  10"2 

9.9581 

0.47316  X  10“° 

0° 

0.0625 

0.71345  x  10"* 

7.1676 

0.65738  x  10“° 

0° 

0.125 

0.91525  x  10"2 

5.5873 

0.84330  x  10~° 

0° 

0.25 

0.13188  x  10- 1 

3.8774 

0.12152  x  lO-* 

0° 

0.375 

0.17224  x  10“ 1 

2.9689 

0.15870  x  lO-* 

0° 

0.5 

0.21260  x  10- 1 

2.4053 

0.19589  x  10-* 

0° 

0.75 

0.29332  x  10" 1 

1.7434 

0.27026  x  10“* 

0° 

1.0 

0.37404  x  10" 1 

1.3672 

0.34464  x  10“* 

0° 

2.0 

0.69691  x  10“ 1 

0.73377 

0.64214  x  10"2 

25° 

-0.1 

0.32506  x  10“* 

27.245 

0.39414  x  10-° 

25° 

-0.075 

0.41572  x  10"* 

21.303 

0.50408  x  10“° 

25° 

-0.05 

0.50437  x  10"2 

17.559 

0.61158  x  KT2 

25° 

-0.025 

0.59348  x  10_if 

14.923 

0.71962  x  10“° 

25° 

0.0 

0.68259  x  10“* 

12.974 

0.82768  x  10“° 

25° 

0.0625 

0.90536  x  10" 

9.7820 

0.10978  x  lO”2 

25° 

0.125 

0.11229  x  lO-1 

7.8870 

0.13616  x  KT2 

25° 

0.25 

0.15650  x  10"1 

5.6588 

0.18977  x  10"  * 

25° 

0.375 

0.20072  x  10- 1 

4.4123 

0.24338  x  lO" 2 

25° 

0.5 

0.24493  x  10" 1 

3.6158 

0.29700  x  lO-* 

25° 

0.75 

0.33336  x  10"1 

2.6566 

0.40422  x  10"* 

25° 

1.0 

0.42179  x  10" 1 

2.0997 

0.51144  x  10-* 

25° 

2.0 

0.77551  x  10“ 1 

1.1420 

0.94034  x  10-* 

50° 

-0.1 

0.45097  x  10"2 

27.149 

0.64296  x  10“° 

50° 

-0.075 

0.54643  x  10"* 

22.406 

0.77904  x  10“° 

50° 

-0.05 

0.64189  x  10"* 

19.074 

0.91514  x  10“° 

50° 

-0.025 

0.73345  x  lO-* 

16.693 

0.10457  x  lO-* 

50° 

0.0 

0.82746  x  10-* 

14.796 

0.11797  x  lO"* 

50° 

0.0625 

0.10625  x  10_i 

11.523 

0.15148  x  lO”* 

50° 

0.125 

0.12975  x  10- 1 

9.4363 

0.18498  x  10-* 

50° 

0.25 

0.17586  x  10_i 

6.9621 

0.25072  x  lO-* 

50° 

0.375 

0.22255  x  10- 1 

5.5014 

0.31728  x  10~* 

50° 

0.5 

0.26922  x  10- 1 

4.5477 

0.38382  x  lO"2 

50° 

0.75 

0.36256  x  10- 1 

3.3769 

0.51690  x  10-* 

50° 

1.0 

0.45590  x  10- 1 

2.6855 

0.64998  x  10“* 

50° 

2.0 

0.82926  x  10- 1 

1.4764 

0.11823  x  lO”1 

72.34° 

-0.1 

0.54127  x  lO-2"" 

28.188 

0.86146  x  10“2 

72.34° 

-0.075 

0.64152  x  10"2 

23.783 

0.10210  x  10“* 

72.34° 

-0.05 

0.74177  x  10~* 

20.569 

0.11805  x  10“2 

72.34° 

-0.025 

0.84202  X  10“2 

18.120 

0.13401  x  10-* 

72.34° 

6.0 

0.94227  x  lO" * 

16.192 

0.14996  x  10“2 

72.34° 

0.0625 

0.11929  x  10"1 

12.790 

0.18985  x  10“* 

72.34° 

0.125 

0.14354  x  10"1 

10.629 

0.22844  x  10-2 

72.34° 

0.25 

0.19309  x  lO- r 

7.9016 

0.30730  x  10“  * 

72.34° 

0.375 

0.24264  x  10 --1 

6.2880 

0.38618  x  10"* 

72.34° 

0.5 

0.29219  x  lO- 1 

5.2217 

0.46504  x  lO" * 

72.34° 

0.75 

0.39129  x  10" 1 

3.8992 

0.62276  x  lO”2 

72.34° 

1.0 

0.49039  x  10" 1 

3.1112 

0.78048  x  lO-* 

72.34° 

2.0 

0.88244  x  10" 1 

1.7290 

0.14044  x  10"1 
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Table  40:  Wetted  perimeter,  mean  velocity,  Poiseuille  number  and  volumetric  flow 
rate  versus  groove  fill  ratio  for  various  values  of  groove  aspect  ratio  (r,*  =  0.0,  0  =  0°). 


p 

w;/w* 

jp* 

V* 

Po 

V* 

0.25 

0.0731191 

8.4006 

0.35142  x  10~8 

0.25 

0.2 

0.223950 

KiMibfeHiHiBl 

9.6365 

0.70676  x  10“6 

0.25 

0.450806 

■igjgfegfeEiiiBai 

0.25 

0.4 

0.734032 

12.059 

0.40386  x  10"4 

0.25 

0.5 

1.04707 

12.514 

0.12238  x  10"3 

0.25 

0.6 

1.36011 

0.29299  x  10"2 

12.849 

0.27336  x  10"3 

0.25 

0.7 

1.64333 

12.871 

0.49958  x  10~3 

0.25 

0.8 

1.87019 

12.940 

0.77004  x  10"3 

0.25 

0.9 

2.02102 

0.25 

0.999 

2.09364 

15.203 

0.18573  x  10-2 

0.5 

0.1 

0.114187 

0.5 

0.2 

0.288528 

8.4052 

0.82718  x  10"6 

0.5 

0.531454 

IHfiM 

8.7511 

0.20062  x  10~4 

0.5 

— 1 

0.827649 

9.3706 

0.13311  x  10“3 

0.5 

0.5 

1.15245 

■rafHSKiaBBi 

9.9581 

0.47316  x  10“3 

0.5 

0.6 

1.47724 

I»ireSEEl»H»Bi 

0.5 

0.7 

1.77344 

0.5 

2.01636 

mwmmtma 

11.282 

0.5 

2.19070 

0.17762  x  10'1 

12.443 

0.5 

0.999 

2.30389 

0.25467  x  10-1 

14.698 

0.12692  x  lO"1 

0.75 

0.1 

0.160025 

0.11137  x  10"5 

10.092 

0.21124  x  10"9 

0.75 

0.2 

0.367435 

0.10959  x  10~3 

8.9207 

0.44516  x  10“6 

0.75 

0.3 

0.634966 

WmWiWSVBM 

8.4299 

0.18655  x  10"4 

0.75 

0.4 

0.951588 

1 

8.5312 

0.17184  x  10"3 

0.75 

1.29509 

8.8579 

0.74738  x  10'3 

0.75 

1.63860 

0.11328  x  lO'1 

9.2444 

0.21238  x  lO"2 

0.75 

0.7 

1.95522 

0.16731  x  10"1 

9.7209 

0.46642  x  10"2 

0.75 

0.8 

2.22275 

ItttklEMRUiM 

■‘TTtntrlT'Ml 

0.75 

0.9 

2.43016 

11.885 

0.17278  x  10"1 

0.75 

0.999 

2.58868 

0.47659  x  lO”1 

13.997 

0.35626  x  10'1 

1.0 

0.1 

0.207697 

MfcblifeiilliM 

10.367 

0.73744  x  10"1() 

1.0 

0.2 

0.453879 

iBBBEa 

9.4362 

1.0 

0.3 

0.752420 

0.66244  x  10~3 

8.7481 

0.13415  x  10"4 

1.0 

0.4 

1.09559 

0.27276  x  10~2 

8.4926 

0.16081  x  10"3 

1.0 

0.5 

1.46369 

8.5590 

0.84758  x  10“3 

1.0 

0.6 

1.83180 

0.12785  x  10"1 

8.8427 

0.27842  x  10"2 

1.0 

0.7 

2.17497 

1.0 

0.8 

2.47351 

0.31143  x  10"1 

10.106 

0.15279  x  10"1 

1.0 

2.71969 

0.46597  x  lO-1 

11.485 

0.32778  x  10"1 

1.0 

0.999 

2.92539 

13.281 

0.69692  x  10"1 
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